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We reformulate the nonperturbative functional renormalization group for the random field Ising 
model in a superfield formalism, extending the supersymmetric description of the critical behavior of 
the system first proposed by Parisi and Sourlas [Phys. Rev. Lett. 43, 744 (1979)]. We show that the 
two crucial ingredients for this extension are the introduction of a weighting factor, which accounts 
for ground-state dominance when multiple metastable states are present, and of multiple copies of 
the original system, which allows one to access the full functional dependence of the cumulants of 
the renormalized disorder and to describe rare events. We then derive exact renormalization group 
equations for the flow of the renormalized cumulants associated with the effective average action. 
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I. INTRODUCTION 

Theoretical investigations of the critical behavior of the 
random field Ising model (RFIM) have a decades-long 
history>ii^ One of their central themes concerns the va- 
lidity of the "dimensional reduction" property, according 
to which the long-distance behavior of the RFIM in d di- 
mensions is identical to that of the system without disor- 
der in d—2. First pointed out by Grinstcinj^ the property 
was shown to emerge from perturbation expansion to all 
orders by Aharony et al^ and Young.— Finally, Parisi and 
Sourlas, in a beautiful 2-page letter)^ related the critical 
behavior of the RFIM to a supersymmetric scalar field 
theory and showed that the supersymmetry leads to di- 
mensional reduction. Phenomenological,i^ numerical^- 
and rigorous studies^— have however established that the 
dimensional-reduction prediction is wrong, at least in low 
dimensions {d< 3). It has also been understood that the 
superfield construction of Ref. loses its validity when 
multiple metastable states are present, which is the case 
in the region of interest.— In addition, it has been shown 
that dimensional reduction follows from supersymmetry 
not only in pcrturbative expansions as in Ref. |6| but in 
a nonperturbative manner f^i"— 

Despite the elegance of the formalism, no steps beyond 
the original superfield formulation have so far proven 
usefuL-ii and the goal of the present work is to provide 
a solution by combining superfield formalism and func- 
tional renormalization group (FRG). This is part of our 
ongoing prograroi^^— to build a consistent and compre- 
hensive theory of the long-distance physics of random 
field models, and more generally of disordered systems, 
based on the nonperturbative functional renormalization 
group (NP-FRG). In papers li^ and ItiS, of this series, 
wc showed that the breakdown of dimensional reduction 
is related to the appearance of a nonanalytic dependence 
of the effective action (or Gibbs free-energy functional 
in the terminology of magnetic systems) in the dimen- 
sionless fields. However, our approach, which was based 
on a "replica method" for handling the average over the 



random field, could not address the pending question of 
supersymmetry and its breaking. (In consequence, there 
was no means to avoid explicit breaking of the super- 
symmetry in the regulators and in the approximation 
scheme.) As the effective Hamiltonian (or bare action) 
of the model in the presence of a random field has always 
multiple minima in the region of interest (near the criti- 
cal point), it would seem that a superfield approach and 
the dimensional-reduction predictions can never be valid. 
Our objective is nonetheless to give a meaning to the su- 
perfield formalism even in cases where multiple minima 
are present and/or supersymmetry is broken, and to in- 
vestigate the validity of the dimensional-reduction results 
as one varies spatial dimension. As will be shown, a res- 
olution of the problem requires an FRG formulation. 

Our starting point is that the fundamental flaw of the 
Parisi-Sourlas supersymmetric construction^ has in fact 
a two-fold origin: 

(1) the presence of metastable states, which can be 
described at zero temperature as the multiplicity of so- 
lutions of the stochastic field equation associated with 
the extremization of the bare action of the RFIM, is not 
counterbalanced by a means to select the ground state 
and 

(2) the use of a single copy of the original disordered 
system does not give access to the full functional depen- 
dence of the cumulants of the renormalized disorder and 
is thereby unable to account for the rare events, such as 
avalanches and droplets, that characterize random-field 
systems. 

The second aspect has already been addressed in our 
previous investigation through the NP-FRG and it re- 
quires introducing copies or replicas of the original sys- 
tem with the same disorder but different applied sources 
(resulting in an explicit breaking of the permutational 
symmetry between replicas). We will show in the com- 
panion paper that the spontaneous breaking of supersym- 
metry (more precisely of "superrotational invariance") 
that comes with the breakdown of dimensional reduc- 
tion is precisely linked to the emergence of a strong 
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enough nonanalytic dependence of the cumulants of the 
renormalized disorder. Curing the first problem on the 
other hand impUcs a way to properly select the ground 
state among all solutions of the stochastic field equa- 
tion. We propose a resolution of this problem through 
the introduction of a weighting factor and the construc- 
tion of a superfield theory in a curved supcrspacc. In 
the present paper, we relate ground-state dominance to 
a formal property that we call "Grassmannian ultralocal- 
ity" . This finding may have value for other problems in 
which a generating functional is built from the solutions 
of a stochastic field equation, as in other disordered or 
glassy systems?^ in turbulence^! or in nonabelian gauge 
field theories.—"— 

The outline of the article is as follows. In Sec. [ll] 
we introduce the superfield formalism for the T = 
equilibrium long-distance properties of the RFIM, as 
proposed by Parisi and Sourlas,— and we discuss the 
symmetries, supersymmetries and the associated Ward- 
Takahashi identities. 

In Sec. mil we discuss the property of "Grassmannian 
ultralocality" and its connection to the fact that a unique 
solution of the stochastic field equation is taken into ac- 
count in the computation of the generating functionals; 
we next formulate a procedure to properly select the 
ground state, which is relevant for the T = equilib- 
rium long-distance properties of the RFIM, through the 
insertion of a weighting factor and the consideration of a 
curved superspace. 

We show in Sec. IIVI that a complete description of 
the renormalized random functional that describes the 
physics of the RFIM, including rare events such as 
avalanches and droplets, requires introducing copies of 
the original system with independently controlled sources 
so that the hierarchy of cumulants with their full func- 
tional dependence can be generated. We then consider 
the properties of the superfield theory with multiple 
copies, its (super)symmetries and the associated Ward- 
Takahashi identitcs. 

In the following section. Sec. |Vl we explore the formal 
consequences of the property of "Grassmannian ultralo- 
cality" . 

Wc next describe in Scc. lVII the generalization to super- 
fields in a curved superspace of our previously developed 
NP-FRG approach in the effective action formalism. We 
carefully discuss the issue of the infrared regulator and wc 
derive the exact FRG equations for the cumulants of the 
renormalized disorder; we also consider the implication 
for these equations of the hypothesis of "Grassmannian 
ultralocality" . 

In Sec. IVIIl wc detail how our formalism allows one to 
describe ground-state dominance at long distance in the 
NP-FRG. Wc show that the limit of infinite curvature, 
which corresponds to a vanishing auxiliary temperature, 
gives back the FRG equations for the cumulants of the 
renormalized disorder under the property of "Grassman- 
nian ultralocality" . 

We finally conclude with a summary and a discussion. 



Additional information is provided in several appendices. 

In the companion paper, we apply the NP-FRG in 
the superfield formalism to describe supersymmetry and 
its spontaneous breaking in the critical behavior of the 
RFIM. We introduce a nonperturbative approximation 
scheme to the exact FRG equations for the cumulants 
and solve the flow equations numerically to determine 
the critical exponents and the fixed-point properties as 
a function of space dimension d. This provides a reso- 
lution of the long-standing puzzles associated with the 
long-distance physics of the RFIM. 

A short account of this work has appeared in Ref . [2^ . 



II. THE PARISI-SOURLAS SUPERFIELD 
FORMALISM FOR THE RFIM 

A. Model and generating functionals 

The starting point is the field-theoretical description 
of the RFIM in terms of a scalar field if{x) in a d- 
dimensional space and a bare action S[ip;h] given by 



(1) 



where = J d'^x, Ub{^) = {t /2)ip^ + {ul 'i\)ip'^ , and h{x) 
is a random source (a random magnetic field in the lan- 
guage of magnetic systems) that is taken with a Gaussian 
distribution characterized by a zero mean and a variance 
h{x)h{y) = As 6^'^\x — y). A (ultra-violet) momentum 
cutoff A, associated with an inverse microscopic length- 
scale such as a lattice spacing, is also implicitly consid- 
ered. 

The Parisi-Sourlas construction goes as followsi^ Tak- 
ing advantage of the fact that at long-distance the ther- 
mal fluctuations are negligible compared to those induced 
by disorder (formally, the critical behavior is controlled 
by a zero-temperature fixed point^ '^^i^'' ), one can focus 
on the ground-state configuration which is solution of the 
stochastic field equation 



5S[v; h] 
5(p{x) 



J{x) 



(2) 



where we have added an external source (a magnetic 
field) J conjugate to the ip field. 

Provided the solution of Eq. ([2|) is unique, the equi- 
librium (Green's) correlation functions of the Lp field are 
obtained from appropriate derivatives of the generating 
functional 



Z,,[J, J] = V^5 



5SbVp\ 
6(p 



-h-J 



det 



X exp / j{x)ip{x), 



(3) 

where Sb[<p] = JJ{l/2)idMx))^ + UBMx))}; the delta 
functional S[ ] enforces that (p is solution of Eq. ([2]) 



3 



and the absolute value of the functional determinant of 
5"^ S b['^]I [5v[x)5ip{y)) is the associated jacobian. Due to 
the postulated uniqueness of the solution, the absolute 
value can be dropped and the functional can be built 
through standard field-theoretical techniques)^ One first 
introduces auxiliary fields: a bosonic "response" field 
(p{x) associated with the integral representation of the 
delta functional and two fermionic "ghost" fields 'ip{x) 
and '4j{x) [satisfying ip'^ = ip'^ — ipip+ipip = 0] to exponen- 
tiate the determinant. We also introduce two fermionic 
sources K{x), K{x) linearly coupled to the ghost fields. 
This leads to 

Zh [J, J, k, K] = exp (Wh [J, J, A% if]) 



TV / VLpV(pVi}Di}cyiY>{ / ~(p{x) 



^ 5Sb[<p\ 
5if{x) 

(pix) [h{x) + J{x)] -f J{x)(p{x) + 'ip{x)K{x) 

where M is an irrelevant constant factor that will be 
dropped in the following and Wh is a random functional 
that depends on the bare quenched disorder (i.e. the 
bare random field h). One then performs the average 
of the partition function (and not of its logarithm as in 
the standard approach to disordered systems) over the 
Gaussian disorder, which provides 



Z[J,J,K,K]^ Zh[J,J,K,K] 

VipVifiVipD-ipcxp { - Sss[(p, <(5, V', + 



(j(x)^(a;) + i^{x)K{x) + K{x){p{x) + J(x)^(a;)) }, 

(5) 



where 



5, 



Vi.x) 



5if{x) 



xjy S(p{x)5ip{y) 



(6) 



We define W[J,J,K,K] = log Z[J, J, K,K], which is 
the generating functional of the connected Green's func- 
tions. An important feature is that due to the identity 
Z[J = 0,J,K = 0,K = 0] = 1, the v3-field connected 
correlation functions of the original problem are obtained 
by functional derivatives of W with respect to J that are 
further evaluated for K = K = J = 0. For instance, 



{ip{x)) {ip{y)) - {ip{x)) {(p{y)) = 



SJ{x)6J{y) 



lK=K=J=0' 



where, since we consider the zero-temperature limit, 
{ifiix)) is equal to the solution of the stochatic field equa- 
tion, Eq. ©. 

The next step of the construction is to introduce a su- 
perspace by adding to the d-dimensional Euclidean space 



with coordinates x = {2;'^} two ant i- commuting Grass- 
mann coordinates 9, 9 (satisfying 9^ = 9^ = 99 + 99 = 
0)1^ By letting x = (x, 9, 9) denote the coordinates, the 
associated (super)metric is given by 

dx^ =dx'^ + ^d9d9 = g,nndx"'dx", (8) 

where {m} = {fi,9,9}, dx^ = dx'^dx^, and a summa- 
tion over repeated indices is implied; the metric tensor 
gmn satisfies: g^^ = (5^^, ggg = ~ggg = -2/ As, with 
all other components equal to zero. With the notations 
dm = d/dx"^,d^ = d/dx^^dg = 8/ 89, etc., one can 
express the corresponding "super-Laplacian" as 



(4) where 



ff""c)„,c>„ = 8^8^ + ABde8g, (9) 
(S™. After introducing the superfield 



$(£) = ^{x) + 9i'{x) + iix)9 + 99ip{x) (10) 
and the supersource 

J{x) = J{x) + 9K{x) + K{x)9 + 99J{x), (11) 
the generating functional can be cast in the form^ 



Z[J] = y 2?$cxp (^-Sssm + y J{x)<^{x)^ 
= exp(W^[^]) 



(12) 



with 



S,M = I [-\<^{x)AssH3l) + UBmx))l (13) 

and with = J^. // d9d9. By a Legendre transform, one 
introduces the effective action r[$] which is the generat- 
ing functional of the 1-particle irreducible (IPI) correla- 
tion functions or proper vertices!^ 



r[$] = -W{J\ + / ^(x)$(ai). 



(14) 



where the (classical) superfield $ and the supersource J 
are related through 



SW\J\ 
5J{x) 



and 



J{x)^ 



5^{xY 



(15) 



(16) 



which can also easily be expressed in terms of the com- 
ponents of the classical superfield, (j), (f>, ip, ?/;, and of the 
supersource, J, J, K. (Note that for simplicity we keep 
the same notation for the superfield $ and its classical 
value, and similarly for the ghost fields we only 

make a distinction for the bosonic fields, i.e. (f> = {ip) 
and (f> = {ip).) 



B. Symmetries, supersymmetries, and 
Ward-Takahashi identities 



The action S'ss[$] is invariant under a large group of 
transformations : 

(i) sign changes associated with a Z2 discrete symme- 
try $ — > — $ (the critical behavior we aim at describing 
is associated with a spontaneous breaking of this Z2 sym- 
metry, i. e. a paramagnetic-to- ferromagnetic transition in 
the language of magnetism); 

(ii) rotations and translations in the d-dimensional Eu- 



clidean space with infinitesimal generators L 



x^d^ and 9^, where ii,v = 1, ...,d, respectively; 

(iii) transformations of the symplectic group with gen- 
erators acting on the Grassmann subspace: i = Ode, 
t = 9dg (associated with "rotations" in the 2-dimensional 
Grassmann subspace) and N = 9d§ — Odg (correspond- 
ing to the "ghost- number conservation"); the generators 
satisfy the commutation relations [t, t] = N, [N, t\ = 2t, 
and [N,t] ^ -2t; 

(iv) translations in the 2-dimensional Grassmann sub- 
space with generators dg,dg: these are linear, BRS-type, 
symmetries)^ 

(v) "superrotations" that preserve the supermetrie and 



can be represented by the generators Qp 



-^ed^ and Q 



The last two sets of transformations mix bosonic and 
fermionic fields and for this reason represent "fermionic 
symmetries" .— They are associated with supersymme- 
tries: the translations in the 2-dimensional Grassmann 
subspace (iv) form a supergroup and the Euclidean ro- 
tations (ii), the symplectic group (iii) and the superro- 
tations (v) form another supergroup known as the or- 
thosymplectic supergroup OSp{2,d)^ 

The linearly realized continuous symmetries and su- 
persymmetries (ii)-(v) lead to a set of Ward-Takahashi 
(WT) identities that can be expressed either at the level 
of the generating functional of the Green's functions or 
that of the effective actioni^ For instance, the invariance 
of the action under the superrotations gives rise to the 
following WT identities: 



J X 



(17) 



and for p > 1 



(18) 

where we have used the notation to indicate a generic 
component of the generator of the superrotation act- 
ing on the superspaee coordinate x; the proper (IPI) ver- 
tices are defined as 



(5pr[$] 



(5$(£i)...(5$(x ) 



(19) 



Similar identites are obtained with the generator (and 
additional ones with the generators 9^. 9^,9^, t, i, 
and N). 

At this point, it is important to stress that it is the 
existence of the superrotation invariance (v) which leads 
to the property of dimensional reduction in the critical 
behavior of the RFIM, with the two Grassmann dimen- 
sions acting as negative dimensions i^iii One should how- 
ever keep in mind that the whole formal construction 
collapses when the stochastic field equation, Eq. ([2]), has 
more than one solution, which is usually the case in the 
region of interest In the following, we show how to 
generalize the Parisi-Sourlas construction in order to pre- 
serve the relevance of the superfield formalism. 



III. "GRASSMANNIAN ULTRALOCALITY" 
AND GROUND-STATE DOMINANCE 

A. On the "ultralocality" of the random 
generating functional in the Grassmann subspace 

Before presenting the method to select the ground- 
state configuration, which should dominate at T = 0, 
we discuss an important property of the random gener- 
ating functional in the superfield framework. To alleviate 
the notations we consider the case of a c? = Euclidean 
space, but the same considerations equally apply to any 
number of Euclidean dimensions. When the stochastic 
field equation [Eq. (2)], which can be rewritten as 



dip 



h+J, 



(20) 



has a unique solution, say ipt{h + J), the generating func- 
tional of the (Green's) correlation functions of the (p field 
is simply expressed as 



Zh{J,J) = e 



J(fi,{h+J) 



(21) 



or, after adding the fermionic sources K, K and using the 
auxiliary fields [see Eq. (4)], as 



Zh{J,J,K,K) = Zh[J] 



Jip,[h^-J)-KKS''^\v,{h+J)y^ 

(22) 

where Sg {(p) is the second derivative of the bare action 
and J" is the supersource introduced in Eq. (|lip . As a 
consequence, the (random) generating functional Wh = 
log Zh is equal to 

Wh[J] ^Jp,{h + J)-KKS^B\ip,{h + J)y\ (23) 

It is straighforward to show that this precisely cor- 
responds to an "ultralocal" form in the Grassmann 
subspace^^ 



Wh[J] = / WhiJii)) 



(24) 
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where 9_ coUects the two Grassmann coordinates 6 and 9 
and Jg = J J dddO; indeed, 



WniJm = JWt'\j) - KKWI^\J) 



(25) 



with Wl^\j) = ip^{h + J) and Wl^'{J) = dip^h + 

J)/dJ = S''s\f*ih + J)y^ as it should be. The whole 
formal construction is fully justified in this case. 

What happens when the stochastic field equation has 
several solutions ? Let us sort the solutions and label 
them with an index a; a = denotes the ground state, 
which is unique for a continuous distribution of the ran- 
dom field, aside from exceptional values of h+J for which 
a coexistence of several ground states may take place. 

The sought-for generating functional that describes 
the equilibrium situation at T = only takes into ac- 
count the contribution of the ground state. As a re- 
sult, one has exactly the same relations as above, with 
ip^,{h+J) replaced by the ground-state solution ipo{h+J). 
Actually, the notation (pQ{h + J) is a little misleading 
as the ground state is generally a piecewise function 
of h + J with discontinuities occuring at special val- 
ues of the latter. These discontinuities are known as 
"avalanches" (or "shocks") and have been observed for 
instance in numerical computations of the ground state 
of the RFIM on 2- and 3-dimensional lattices,^"-^^ The 
ground-state solution ipQ{h + J) should rather be written 
as J2iCo,i{h + J)'fns{h + J) with Co, i the characteristic of 
the ith interval along the axis h + J in which the ground 
state ifo.i is a continuous function oi h + J (see also Ap- 
pendix |A| . More generally, including a single solution in 
the generating functional Zji{J, J, K, K), provided this 
solution is piecewise defined for all values oi h + J, also 
leads to the "ultralocal property" in the Grassmann sub- 
space for yVh[J], Eq. (p4|. 

On the other hand, one easily realizes that this prop- 
erty cannot be true when several different solutions are 
included in the generating functional, no matter how one 
chooses to weigh these solutions. The weighting of the 
solutions that corresponds to the above Parisi-Sourlas su- 
persymmetric construction, Eq. is given by the sign 

of the determinant of the Hessian Sg'^i^ Then, setting 
to zero the fermionic sources for simplicity, one has 

Z„(J, J) = ^(-l)"(")C„,,(/i -t- J)eJ^''Ah+J)^ (26) 



(2), 



where n(a) is the index of the ath solution, i.e. the 
associated number of negative eigenvalues of the Hessian 

(2) 

Sg , and Ca,i{h + J) is the characteristic function of the 
ith interval over which ipa = ^Pa,i is a continuous function 
of ft, -I- J [at the boundaries of Ca,iih + J) the solution 
(fa has a discontinous jump or possibly stops to exist]. 
Obviously, 



W„(J, J) =log^(- 



(27) 



where (pa denotes J2i^a,iih + J)<PoL,i{h + J), is not "ul- 
tralocal" in the Grassmann coordinates, i.e. cannot be 
put in the form of Eqs. P3II25|) [observe in particular that 
the form in Eq. (j23p implies a linear dependence on J, 
which is not satisfied by the above expression]. When 
J = 0, one finds that Zh = 1 and Wh = 0, just as in 
the case where a unique solution is taken into account)^ 
However, the result now follows from the property that 
the sum of all solutions weighted by the sign of the de- 
terminant of the Hessian is a topological invariant which 
is equal to 1 in the present casei^ 

Note that contrary to the toy model discussed in 
Ref. |34| , in which the breakdown of the supersymmetric 



formalism is related to the fact that the stochastic equa- 
tion has no solutions for a certain range of the random 
field, the stochastic equation associated with the present 
scalar field theory in a Gaussian distributed random field 
has always at least one solution (the ground state always 
exists). 



B. Selecting the ground state 

A natural procedure to select the ground state is to add 
in the generating functional, Eq.Q, a weighting factor 
that strongly favors the solution with the smallest ac- 
tion. This can be done through a Bolzmann-like factor, 
namely. 



Vif 5 



6ip 



h- J 



dot 



StpStp 



exp 



13{Sb[^]- I [J{x)+h{x)Mx)) + J J{xMx)] 

(28) 



where /3 is the inverse of an auxiliary temperature (the 

actual temperature is equal to zero). Note that zjf'' is 
not the same as the partition function obtained from the 
equilibrium Boltzmann-Gibbs measure at a temperature 
T = Even if /? is large enough that only minima 

contribute to Eq. (PSI) , the latter expression only includes 
the contribution of the minima whereas the Boltzmann- 
Gibbs measure also takes into account the contribution of 
the basins of attraction of the minima (roughly speaking, 
the thermal fiuctuations around the minima). 

With the above generating functional, one finds that 
the average of the field ip is given by 



JlogZ;;' 



SJ{x) 



1 



J=0 



[J = 0, J] 5J{x) 



J=0 

(29) 

The simplicity of the Parisi-Sourlas formalism is however 
lost as zjf\j = 0, J] 7^ 1. In consequence, simply con- 
sidering the average over the disorder of zjf^ is no longer 
sufficient to generate the (/3-field correlation functions of 
the original problem. 

Interestingly, the superfield formalim can still prove 
useful. After having introduced one bosonic and two 
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fermionic auxiliary fields as before, added two fermionic 
sources, grouped all the sources in a supersource as in 
Eq. (jlip. and similarly grouped all the fields in a su- 
perfield as in Eq. ([TU|) . the generating functional can be 
rewritten as 



dede[i + i3ee]SB[^{0)\ 



dede[i + i3ee] [h{x) + j{x,e)] 



= exp(wf)[j7]). 

(30) 

The construction may appear rather formal and un- 
tractable but it turns out that Eq. (|30l) can be expressed 
in a way which will prove efficient to study the symme- 
tries of the theory and investigate its long-distance prop- 
erties. It is convenient to introduce a superspace combin- 
ing the d Euclidean and the 2 Grassmannian dimensions 
in which the Grassmann subspace is now curved. To be 
more specific, we replace the metric tensor of the Parisi- 
Sourlas formalism (see section II-A) by 



(1 



909 — See — : 



(31) 



keeping gmn ~ Sp,i, for the Euclidean sector and all cross- 
components between Euclidean and Grassmannian coor- 
dinates equal to zero. So long as we are not interested in 
mixing bosonic and fermionic directions (mixing occurs 
when considering superrotations, see section II-A), there 
is no need to introduce the factor 2/ As in the superme- 
tric as was done in the Parisi-Sourlas formalism above. 

The properties of the curved superspace are discussed 
in Ref. |37| . In a nutshell, one has the usual prescrip- 
tions of Riemannian geometry that in order to satisfy the 
isometrics of the curved Grassmann subspace, one should 
(i) contract Grassmann indices with either the metric 
tensor or its inverse, (ii) integrate over the Grassmann co- 
ordinates with a measure y^sdetg dOdO = (1 + P09) d9d9, 
where sdetg = (1 -I- (399)'^ is the superdeterminant of the 
metric tensor in the Grassmann sector, (iii) use, if neces- 
sary, covariant derivatives as well as the proper Laplacian 
operator Ae_ — .g™" [dmdn ~ T^^dp), where m,n,p = 9, 9, 
summation over repeated indices is implied, and the 
r^„'s are the ChristofFel symbols with nonzero compo- 
nents = -r% = (39 and T^g = -T^ = 139; the 
Laplacian for the Grassmann subspace is then explicitly 
given by 

Ae = 2(1 -I- I399)idedg - 0dg - p9de). (32) 

It is easy to show that the parameter f3 (inverse of an 
auxiliary temperature) is up to a factor 1/6 equal to the 
Ricci scalar curvature of the Grassmann subspace 

We now come back to the discussion of the previous 
subsection concerning "ultralocality" in the Grassmann 
subspace. For ease of notation, we consider again a d = 
Euclidean subspace and we momentarily drop the two 



fermionic sources. The random generating functional 
y^jf"^ ['J] corresponding to the superfield construction in 



Eqs. (f29| and (|30|) is expressed as 



wr(j,j) = 

- l3[SB[Mh + J)] - (h + J)Mh + J)] + JMh + J) + 
log (l + ^(_i)"(")e-^'^'S'„,o(/i+J)+J[vo(/»+J)-yo(/»+./)]^^ 

(33) 

where ASafl = [SsVPa] - Sb[^o]) - [h + J){(pa - ^o)- 

Assume first that the parameter /? can be taken suffi- 
ciently large for ensuring ASa,o for all solutions 
a 7^ 0, so that all contributions but that of the ground 
state can be neglected in Eq. l\'S'3\i . It is then easy to show 
that the random functional can be put in an "ultralocal" 
form appropriate for the curved Grassmann subspace. 



(34) 



where 



M^f ^ [J] - Sb [ifoih + J)] -{h + J)ipo{h + J) (35) 

is actually independent of (3 and where the integral over 
9_ now involves the metric factor (1 + 1399), i.e. Jg = 

J J{1+ I399)d9d9. 

In fact, no matter how large (but finite) /?, the above 
assumption may not be valid for all realizations of the 
random field. Excited states with ASa,o ^ 1//? may 
be present. The vast majority of them involve only lo- 
cal changes of configuration with respect to the ground 
state.— 1^ As a result, the term J{ipa — fo) remains 
small in d Euclidean dimensions, provided J of course 
stays small or finite (here, we are ultimately interested 
in the limit J = 0). In addition, rare excitations, also 
known as "droplets",— may involve a large-scale reor- 
ganization of the ground-state configuration while satis- 
fying AS'q^o ^ and therefore provide a significant 
contribution to wjf'^. In such a situation, the "ultralo- 
cality" in the curved Grassmann subspace is broken and 
the random generating function has the general form 



Wt'\ji9),il 



:99)de_J{9),AgJ{9)], 



(36) 



where de is a short-hand notation for designating either 
dg or d0~a.nd we recall that Ag_ = 2{1+I399){dgdg-I39dg- 
(39 dg) is the Laplacian in the curved Grassmann subspace 
(see above). Note that the functional dependence on the 
Euclidean coordinates is completely general at this point. 

The droplets being nonetheless rare events, and the 
other excitations being essentially local, we expect that 
the deviation from "ultralocality" is small when j3 is 
large, i.e., 



W, 



Wlf''^^[j{0)] + corrections. 



(37) 
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with the corrections going to zero (and wjf^'^^ going to 
a well defined limit) when /3 — >■ oo. In the following, we 
will check the correctness of this behavior and show how 
in the FRG flow the contributions coming from errors in 
selecting the ground state actually become subdominant 
as one approaches the critical fixed point. 



IV. CUMULANTS OF THE RENORMALIZED 
DISORDER AND THE NEED FOR MULTIPLE 
COPIES 



the system is considered (see section |TT| and in the latter 
the sources acting on the replicas are all taken equal. As 
a consequence, in both cases, one has only access to cu- 
mulants in the specific configuration with all arguments 
equal. Quite differently in the present formalism, we con- 
sider multiple copies or replicas and supersources that ex- 
plicitly break the (permutational) symmetry among these 
rcplicaS)^ We note in passing that, by construction, this 
takes care of the problem coming from having to aver- 
age the logarithm of the partition function over disorder 
when /3 ^ 0. 



A. Why the need for multiple copies ? 

A central quantity is the random ( "free energy" ) func- 
tional yV^'^''[J'], introduced above. This random func- 
tional is characterized by its (functional) probability dis- 
tribution or, alternatively, by the infinite set of its cu- 
mulants (if of course the cumulants exist). Dealing with 
cumulants has the advantage of involving an average over 
the bare disorder: as a result, one recovers the trans- 
lational (and rotational) invariance in Euclidean space 
which is otherwise broken by the space-dependent ran- 
dom field. In the following, we will therefore consider a 
formalism based on cumulants. However, a crucial point 
when working with such disorder-averaged quantities is 
that one does not want to lose track of the rare events 
that characterize systems with quenched disorder. For 
random-field models, these rare events are expected to 
take the form of "avalanches" or "shocks" that are seen in 
the dependence of the ground state on the applied source 
at zero temperature and of low-energy excitations known 
as "droplets" at nonzero temperature (see above). As 
shown in previous workfi^^— i^"— these phenomena show 
up in the cumulants of the renormalized disorder as a sin- 
gular dependence on the arguments (for an illustration in 
a simple zero-dimensional toy model, see Appendix |^ . 
Describing such features therefore requires the functional 
dependence of the cumulants for generic arguments. For 
instance, a complete characterization of the random func- 
tional Wjf'' [JT] implies the knowledge of all its cumulants. 



defined as 



;('3)r 



which are 



(38) 



B. Multi-copy superfield formalism 

The cumulants of yvjf^ for generic arguments can be 
obtained from the average over the bare disorder of the 
extension of Eq. (|30p to an arbitrary large number n 
of copies submitted to independently controllable (su- 
per) sources, namely. 



exp(w(^)[{J-4]) = nz(«[J-J 



exp(5^wf)[J.]). 



(40) 



With the help of the curved superspace introduced above 
and by combining Eq. ([50]) and Eq. PO)) . we end up with a 
superfield theory associated with the following partition 
function: 

2(^)[{Ja}]-CXp(W(^)[{J-4]) 

= / np'&aexpf'-^tW [{$,}] +f] / J'„(x)$,(£)), 

where the multicopy action is given by 
" r 1 



(41) 



71 n 



^ilT. I I 'fa,(x,^i)$.,(x,^2) 



(42) 



Wf'[Ji,J2] 



(39) 



etc. Generic^ i.e. independently tunable, arguments re- 
quire the introduction of several copies or replicas of the 
original system, each with the same bare disorder (ran- 
dom field) but coupled to different external supersources. 
It is worth stressing that this is not what is done in the 
Parisi-Sourlas supersymmetric approach nor in the con- 
ventional replica trick. In the former, a single copy of 



and it should be kept in mind that the integral over the 
Grassmann coordinates include a metric factor due to 
the curvature /3 (as a resuh, = / + I360)ddd0). 
Note that for n >2 the first (kinetic) and last (disorder- 
induced) terms in the above expression of the multicopy 
action can no longer be combined and simply expressed 
with the super-Laplacian, as it was the case for the 1- 
copy action S'55[$] in the absence of curvature (/3 0) 
which is defined in Eq. ()13|) . 



The cumulants of W^^^ can now be generated by the 
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expansion in increasing number of sums over copies, 

n n ^ 



p>l ai — l dp — 1 

(43) 

where the pth cumulant Wp^' is fully symmetric under 
any permutation of its arguments (and independent of 
n). As is obvious from the above equation, at least p 
distinct copies must be considered to describe the pth 
cumulant with generic arguments, so that an arbitrary 
large number of copies (or replicas) are needed to gener- 
ate all cumulants. 

By using Eq. p6p. the first cumulant can be formally 
rewritten as 



(44) 



where W^^\l] is a short-hand notation for indicat- 
ing a functional of Ji(^i), (1 + f ^16*1)% Ji(^i), and 
iTi (^1 ) (the functional character comes from the Eu- 
clidean dependence that remains completely general at 
this point); there is no additional explicit dependence on 

The second cumulants 



Grassmann coordinates in W{ 
reads 



(45) 



with an obvious extension of the above notation, and sim- 
ilar expressions hold for the higher-order cumulants. A 
physical interpretation of the cumulants when the ran- 
dom functional W^'^'' is "ultralocal" in the Grassmann 
coordinates will be given in the next section. 



C. Legendre transform and effective action in a 
curved superspace 

Due to the specific form of the source term in the 
presence of curvature, i.e., explicitly, Ix 1 1 d9dd{l + 
pee)^a{x,9)Jaix,e), the functional W^'^'>[{Ja}] does not 
exactly generate the average of the superfields; rather, 
one has an extra metric factor, 

.V = i^+m^a{x,e). (46) 
0ja{x,9) 

The Legendre transform defining the effective action is 
then expressed as 

(47) 

From this equation, one can determine the relation be- 
tween the second functional derivatives F^^' and W'^) 
(where we have omitted the superscript (/3) to alleviate 



the notation, as we shall do each time superscripts indi- 
cating functional derivatives are involved): 



5^x2,12) 



= 5ab5^'^\xi - X2)59^e^ 



(48) 



so that 



{1 + /39iei)SabS<-''Hxi - X2)Sg^e^ 



E 



(1 - 2^0303) 



9„ Jx3 



5Ja{xiAi)5Jc{x3,e^) 5^,{x3,e_^)5<^b{x2A2y 

(49) 

where 5b_^i^ ^ Sg^gje^e^ = (6*1 - 92){di - 6*2) (with this 
definition, Sg^g^ = l + pOiBi). F^^) and W^^) are thus 
essentially inverse operators, provided that the effect of 
the curvature of the Grassmann subspace is appropriately 
taken into account. 

Like the generating functional W*-'^-' [{J7a}], the effec- 
tive action F'^''^ [{<I>a}] can be expanded in increasing 
number of unrestricted sums over copies (considering now 
the superfields {^a} SiS fundamental variables in place of 
the supersources {Ja})' 

r(«[{<i>a}] = E E - E 



p>l ai—1 ap — 1 



p\ 



(50) 

where V^'^ is a fully symmetric functional of its p argu- 
ments whose functional form is independent of the num- 
ber n of copies. The sign (— 1)^^^ is chosen for further 
convenience. 

The above expansion is similar to that in number of 
unrestricted (free) replica sums developed in Ref. [l8| . 
In consequence, one can apply, mutatis mutandis, the 
results of Ref. [l^ and relate the v'^p^ 's to the cumulants 
Hp^'* as follows [we drop for simplicity the superscript 
(/3) in the expressions]. The first-order term r^''''[$] is 
the Legendre transform of yV^'^''[J']; namely. 



with 



(51) 



(1 + i3ee)^{x) = ^J^iE^ = (52) 



8J{x) 

whereas the second-order term r2'^^[$i,$2] is given by 



(53) 
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where is the nonrandom source defined via the in- 

verse of the Legendre transform relation in Eq. (|52p . i.e., 



(54) 



The above expression of V2 motivates our choice of signs 
for the terms of the expansion of the effective action F, 
Eq. ([50|) : r2[$i,$2] is directly the second cumulant of 
'WhiJ] (with the proper choice of J[^]). 

For the higher-order terms, one finds after some 
lenghty but straightforward manipulations (see Ap- 
pendix |B| 



(2) 



X W['°' j^[$i],^[$3]] +per?7i(123) 



(55) 



where perm(123) denotes the two additional terms 
obtained by circular permutations of the superfields 
$i,$2,$3 and V\ ' is the inverse of W-j; ' through an 
equation similar to Eq. (P^. and so on. This procedure 

leads to a unique functional form for fp^^ which is ex- 
pressed in terms of cumulants of "W^^ of order y or less. 
(Note that this guarantees that the functional form is 
independent of the number n of copies.) 

As illustrated by Eq. T^''^ $p] for p > 

3 cannot be directly taken as the pth cumulant of a 
physically accessible random functional, in particular 
not of the disorder-dependent Legendre transform of 
W,j^''[J'[$]]. However, as it can be expressed in terms of 
such cumulants of order equal to or lower than p, we will 

call for simplicity the fp'^^ 's "cumulants of the renormal- 
ized disorder" (which is true for p = 2) in what follows. 
See also section V-B below. 



D. Symmetries and WT identities for multiple 
copies 

The presence of curvature in the Grassmann subspace 
and of multiple copies still allows invariance of the theory 
under a large group of transformations. The multi-copy 
action, which is given in Eq. (|42p . is indeed invariant 
under the following symmetries: 

(i) The permutations 5„ among copies and a global Z2 
symmetry; 

(ii) the global rotations and translations in the d- 
dimensional Euclidean space; 

(iii) the symplectic transformations with generators 
t = Qde , t = ddg and iV = Qdf) — Qde acting on the 2- 
dimensional curved Grassmann subspace, independently 
for each copy; 



(iv) the two isometrics of the curved Grassmann sub- 
space that generalize the translations of flat space, inde- 
pendently for each copy; their generators are (1 — fiBQ)de 
and (1 - ped)dg. 

(v) It was shown in Ref. [s^] that the above isometrics 
(iii) and (iv) are the only possible ones in the presence 
of a nonzero curvature /3. However, when the curvature 
is set to zero and, in addition, when restricting the su- 
persources such that in all copies except a given copy a, 
the components Ji, = Kb = Kf, = {b a), the partition 
function in Eq. ()41|) is invariant under the superrotations 
considered in section HFB. In this case, one is effectively 
back to a 1-copy system since it is found that 



J V^a CXp [ - Sssl-i-a] + Ja{x)<^ a{x)\ , 



(56) 



where Sss is given in Eq. (|13p and where we have used 
that the 1-copy partition function Z'^''=°)[j76 = Jh, h] ~ 1 
for all copies b ^ a. Invariance under the superrotations 
follows directly. 

The above continuous (super) symmetries imply a set 
of WT identities satisfied by the generating functionals 
W(^)[{Ja}] and F(^)[{$a}]- Denoting by Vb_ any one of 
the generators acting on the Grassmann coordinates 9_ in 
a chosen copy a, one finds that 



W^^^[U/}])=0, (57) 



where we have again dropped the superscript (/3) to 
avoid confusion with the superscripts indicating func- 
tional derivation. This WT identity carries over to the 
effective action: 



1 



T^}[{^f}]VeJ>a{x) = Q, (58) 



where we recall that the integral over the Grassmann 
coordinates comes with a factor (1 + pOO). 

Through differentiation with respect to the superficld 
$a, the above equation leads to WT identities for the IPI 
vertices. For the symplectic transformations, by using 
the fact that Eq. (|58p can be rewritten as 



(1 

one finds for p > 1 
p 



$a(x)fFW [{$/}] =0, (59) 



9=1 



+ / $a(£p+i)Wir[Cii...K,,.^^^) [{$/}] = 0, 

where tq = Oqdg , and similarly for the transformations t 
and N. 
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For the generalized translations (1 — I369)de and (1 — 
/399)dg, a little more care is needed; Eq. ([55]) can now be 
reexpressed as 



1 



'i>aix)de{[l - /30"e]rW [{<&/}]) = 0, (61) 



so that the WT identities for the IPI vertices become 



fl 



l\p(p) 



9=1 



)...{apX ) 



$a(£,+i)a,,,ji-/?0,+i0p+i]r[^+\\ (^^^^^ ^ J [{$/}] = 

.(62) 

and similarly with the other generator. These WT iden- 
tities generalize those already encountered at the 1-copy 
level for a flat Grassmann subspace in Sec. II-B. 

In addition, when the curvature (3 is set to zero and the 
supersources are restricted as discussed above, the super- 
rotational invariance of the partition function also leads 
to WT identities. After making use of the invariance un- 
der translations and symplectic transformations in the 
(flat) Grassmann subspace, one finds that the solution of 
the Legendre relation for a supersource J7h(x) = Jb{x), 
i.e.. 



[{*/}] 



(63) 



where {$/}' denotes the set of all copy superfields but 
satisfies $;,(x) = 4>b{x) (with ipb{x) = ipb{x) = (/>f,(a;) = 
0); one then derives the following WT identities for the 
superrotation invariance: 



(64) 



where the field components ipb{x), ipbix), (j)bix) have been 
set to zero in all copies b ^ a and Qx is defined in sec- 
tion II-B. A similar expression holds with the generators 
Qx- By functional differentiation, WT identites are also 
obtained for the higher-order IPI vertices. This will be 
further exploited in the next section. 



V. EXPLORING THE FORMAL 
CONSEQUENCES OF "GRASSMANNIAN 
ULTRALOCALITY" 

A. Expansion in "ultralocal" cumulants 

Assume now that the random functional wlf"^ is 
"ultra-local" in the Grassmann coordinates, which, as 
discussed before, implies that a unique configuration is 
included in the computation of the random generating 
functional for each realization of the random field. In 
this case, the expansion of the generating functional 
W'^^[{j7a}] in increasing number of sums over copies. 



Eq. p3|) . coincides with a "multilocal" expansion in 
Grassmann coordinates, i.e., 



W^^'^[Ji,...,Jp] = 



w^^^[Me,),...,Me^)], 

(65) 

where Wp^^ no longer depends on the derivatives of the 
supersources in the Grassmann directions and is the pth 
cumulant of the "ultralocal" functional W^^\j'{9)] de- 
fined in Eq. (p4)) . Note that due to the assumed unique- 
ness of the solution included in the computation of the 
random generating functional, M^ft,[J'(^)] as well as its 
cumulants, 



Wi[Ji{9,)] = WH[Jii9,)], 



W2 [Jl {9, ) , ^2 (^2 )] =Wh [Jl ie,)]Wh [J2 (02)] 



(66) 



(67) 



etc, are independent of the curvature /3 (see section III- 
B). By an abuse of language, we will characterize the 
cumulants obtained from an "ultralocal" random func- 
tional as "ultralocal" in the following. 

A physical interpretation of wjf'' and its cumulants 
is next obtained by restricting the supersources to their 
physical component, J{x) = J{x), which plays the 
role of an applied magnetic field. [Said otherwise, this 
amounts to considering supersources that are uniform 
in the Grassmann subspace, i.e. to set 8 = 9 ~ in 
the defining expression in Eq. (|lip .] Wh is then given 
by Eq. ([55]) (properly generalized to d-dimensional Eu- 
clidean space) and, as stated above, is independent of the 
auxiliary parameter /3. Its first derivative is by construc- 
tion equal to the ground-state configuration </Jo[/i + J]. 
The first cumulant Wi[J] then gives access to the ther- 
modynamics and its first derivative is the average of the 
physical field (the "magnetization"), which corresponds 
at zero temperature to the ground-state configuration. 
Its higher-order derivatives 



6PWi[J] 



5J{xi)...5J{xp) 



(68) 



correspond to Green's functions that are related to what 
is known in the literature on disordered systems as "con- 
nected" correlation functions of the field. For instance, 
the second derivative is related to the linear response of 
the magnetization to a change of the applied magnetic 
field. 

The higher-order cumulants describe the distribution 
of the renormalized disorder. They generate through dif- 
ferentiation, and after setting all sources equal (Ji — 
J2 = • • • = J), the so-called "disconnected" correlation 
functions of the original system. We introduce the nota- 
tion 



51^+-+1.Wp[Ji...Jp] 



(69) 



5 J i{x 11).. 5 J i{xiq^)... 5 J p{x pi).. 5 J p{Xpq^) 
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Then, for instance, W2^]:^ [J, J] is equal to the standard with 
two-point "disconnected" correlation function defined in 
Eq. ([T]), which, in magnetic systems, is directly accessible 
to experimental measurementsi^ 



B. Expansion of the effective action 

From the above equations and the definition of the 
Legendrc transform, one derives that the effective action 
has also a "multilocal" expansion in Grassmann coordi- 
nates, which corresponds to having an "ultralocal" form 
for each term of the expansion of F*^^^ in number of copies 
[Eq. dSni)], i.e., 

<i>,J= / ... / T<f^M6i),-,mp)], 
J0, Je 



(70) 

where Tp^'^ does not contain any additional explicit de- 
pendence on Grassmann coordinates (nor any depen- 
dence on the derivatives of the superfields in the Grass- 
mann directions). As before, the dependence on the Eu- 
clidean coordinates, which actually still makes the Fp'^'''s 
nonlocal functionals of the superfields, is left implicit. 
The proof is easily derived from the expression of the 
rlf'^'s in terms the Wq^'^'s with q < p and the "ultralo- 
cal" property of the latter. A more specific discussion of 
the relation between the Tp^'^ 's and the wj/'^ 's is provided 
below. 

From the above property one can derive expressions 
for the proper (IPI) vertices, defined as 



r 



^Pr(^)[{$4] 



(71) 



which will prove useful in the following. After specializing 
to fields in the physical subspace, ^a{x) = (f'aix), which 
are relevant for the equilibrium behavior of the RFIM, 
and introducing the notation 



^ p-xii..xig-^ 



[01. 



j9i+...+gpr„ 



H...9pl 



(72) 



S(j)i{xn)-.S4>i {xiqj...6(f>p{xp^)..6(j)p{xpqj' 



one obtains the following expression for the 1-point 
proper vertex. 



rSaU)[{<^'^}] - (1 + /3^iei)|ri;ij0aJ - /3 X 

and for the 2-point proper vertex, 



(73) 



02 

where Sg_^0^ is defined below Eq. (|49)) . and 



(75) 



2;xi,X2 ' ^0-2] 



(76) 

As the order increases, the formulas go along the same 
lines but become more involved; for instance, one finds 

r(aiXi),(a2X2).(a3X3)[{<^''}] = '^"10203(1 + PQlOl)50^e^9^ 



l;xiX2X3l'f^aiJ / 2:xia 

5„,„,(l+/35l0l)(l + /3^202)<5fl^e,|rgj,2,,3[<^ai,0a 

E r3;xl'x2,x3,. [<^ai > '^a3 , 0aj + • • • j> + perm(123) 

(1 + /3^i0i)(l + fd^d^)(\ + fOzOz)\^t^,,,^xM-^ > 'l>a,Aa,\ 

(77) 

etc, where (Je^e^e, = ^l^lj'l^l^ i^^d perm( 123) denotes 
the two additional terms obtained by circular permuta- 
tions of the indices 1,2,3. 



C. Interpretation of the Fp's 

The relations between the Fp's and the Wp's follow 
from Eqs. (|51ll55p and Eq. They are straighforward 

for the first terms, but get more involved as the order 
increases. To obtain full information, it is sufficient to 
consider configurations of the superfields that are uni- 
form in the Grassmann subspace. (Note that the Wp's 
being independent of /3, so are the Fp's.) 

More precisely, one obtains that Fi[(/)] is the Legendre 
transform of W\ [J] , namely, 



Fi[0] = -VFi[J]+ / J{x)^{x) 



(74) 



with 

The second-order term is given by 

F2[(/)i,02]=VF2[J[</'l],J[</'2]], 



(78) 

(79) 
(80) 



where J[(j)] is the (physical) nonrandom source defined 
via the inverse of the Legendre transform relation in 
Eq. (HSl), i.e., 

j[0](x) = r«[0]. (81) 

and the third-order one by 

r3[01,<^2,03] = ~W3[J[(f>l],J[(f>2lJm] + 

I (<^.'[J[</<i],J[</<2]](M^f'[J[</.i]])"' 

J xy L ^ xy 

X M-yl[J[<Pi].J[H] +perm(123)|, 

(82) 

etc., where perm(123) denotes the two additional terms 
obtained by circular permutations of the fields 0i, 02, 03- 
As stated above, by an abuse of language, we will gener- 
ically call the Fp's "cumulants of the renormalized disor- 
der". 

As we did in Ref. |18| , it is also instructive to introduce 
a "renormalized random field" /«[(/)] (x) as the derivative 
of a random free-energy functional (it can equivalently 
be defined at the level of superfields), 

mix) = {wu[jm - Wh[m]) , (83) 

where J[(/)] is the nonrandom source given by Eq. (|8ip . 
The first moment of h[(f)\ is equal to zero by construction, 
and it is easy to derive that the pth cumulant [p > 2) 
is given by the derivative with respect to (/>!,..., 0p of 
Wp[J[0i], this derivative in turn can be related 
to derivatives of the F^'s with q < p; for instance. 



M<^i](xi)/i[(/)2](x2) = r(;,^|..J0i,02]. (84) 

The cumulants of order p > 3 are given by 

^jy,xu--,xp[4'ii (^pl plus additional terms involving 
higher derivatives of F^'s with q < p; again, by an abuse 

of language, we will simply refer to T^'"^^ as the "pth 
cumulant of the renormalized random field" 

D. WT identities for the superrotational invariance 

Finally, we make use of the above developments to de- 
rive explicitly the WT identities associated with superro- 
tational invariance when /3 — and the multi-copy the- 
ory is reduced to a one-copy problem by an appropriate 
choice of the supersources (see above). We start with 
Eq. ([M)) . which we functionally differentiate to introduce 
higher-order IPI vertices, and we consider configurations 
of the superfields that are uniform in the Grassmann sub- 
space but nonuniform in the Euclidean subspace. Thanks 
to the expressions in section V-B, we decompose each 
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identity into components associated with different poly- 
nomials in the Grassmann coordinates. We then find sev- 
eral types of relations: first, relations merely expressing 
the translational invariance in Euclidean space, namely, 

diA%[<l>] = -f '^(^2)52^r(^l,,J<A], (85) 

•1x2 

Jx3 

(86) 

(9i^ + 92p)F(;,^|,,J0,0] = 

Jx3 

(87) 

etc. 

Secondly, we also obtain more specific and interesting 
identities that relate Fp and Fp+i, e.g., 

di^^Zx. ^] -^K- 4)ri;i,., M = 

r (88) 

- / 0fe)93^rg|,3_,j</,,0], 

J X3 

5iMr3"iL,:.3[?i»] - 

(89) 

etc. For a uniform physical field 4>{x) = (j), Eq. ([55)1 
becomes 

di.^^^lxA^^^) - ^(^1 - = 0, (90) 

which after Fourier transforming and using the transla- 
tional and rotational invariance in Euclidean space leads 
to 

F(")(g2;^,^)=ABa,.r(2)(g2;^), (91) 

with the obvious notation: F^"^^ ,^^ (<?!), (/)) = {27r)'^ 5'^'^\qi + 

q2)T^^\q1;4',(l)), etc. Similar identities are derived for 
the higher orders. 

VI. NP-FRG IN THE SUPERFIELD 
FORMALISM 

A. Nonperturbative FRG 

The main purpose of the present work is to develop 
a formalism that allows one to study the supersymme- 
try, more specifically the invariance under superrotations. 
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and its spontaneous breaking in the RFIM. To this end 
we upgrade our NP-FRG approacbi^"— to a superfield 
formulation and use the tools developed in the previous 
sections to extend the Parisi-Sourlas formalism to the 
relevant situations in which multiple minima of the bare 
action may be present. The key points involve: 

(1) adding an infrared regulator that enforces a pro- 
gressive account of the fluctuations while ensuring that 
the initial condition of the RG flow satisfies the (super) 
symmetries of the action in Eq. and corresponds to 
a unique solution of the stochastic field equation, 

(2) considering copies of the original disordered system 
that give access to the full functional field dependence of 
the renormalized cumulants of the disorder, 

(3) selecting the ground state through the introduction 
of a proper weighting factor and use of a curved super- 
space, 

(4) using the WT identities associated with the (super) 
symmetries to ensure that neither the regulator nor the 
approximations explicitly break the latter. 

Extending our previous work to the superfield theory, 
we introduce a generating functional of the correlation 
functions at the running scale k for an arbitrary number 
n of copies of the system (coupled to the same random 
field but submitted to different external sources) and a 
weighting factor involving the auxiliary parameter /3, 

« n 

4'\{Ja}]^ J n^'I'aexpl-A^f [{$4] 

(92) 

-5(«[{$4] + E / M^D-faix))}, 

where S*^^^ [{^a}] is defined in Eq. (|42|) and involves 
the curved superspace measure as in the preceding sec- 
tions. As previously discussed, the n-copy action is in- 
variant under the Sn permutational symmetry, the Z2 
symmetry, the translations and rotations in the Eu- 
clidean space, and, separately for each copy, under the 
isometrics of the curved Grassmann subspace. The reg- 
ulator is as usual taken quadratic in the superfields. De- 
manding that it satisfies the above symmetries and re- 
quiring in addition that it keeps the multilocal form of 
the bare action in Eq. (j42[) imply the following form: 

^^r = j 11 I I *«ifel)7^Mia2fe,^2)^a2(£2): 
a\ .a2 — -L — 1 — ^ 

where TZk.aia2 denotes infrared cutoff functions satisfying 

T^k,aia2{xi,x2) =<5aia2<5eje2 (1 " l3di9i)Rki\xi - 0:21) 
-f- Rki\xi ~ X2I), 

(94) 

where Se^g^ is defined below Eq. P^ . The infrared cutoff 
functions are chosen such that the integration over modes 
with momentum |g| <C is suppressed (see below); these 
functions must go to zero when fc — >■ so that full inte- 
gration is recovered in this limit J^ii^iS [More precisely. 



we shall see below that Rk{q^) goes to zero except for 
(7 = 0, which nonetheless does not alter the property 
that the regularized theory converges to the full theory 
when k — > 0.] 

The central quantity of our NP-FRG approach is the 
so-called "effective average action" F^.'^'' which is the 
generating functional of the IPI (proper) vertices^^ at 
scale k and is obtained from W^^'' = log z\^^ by a mod- 
ified Legendre transform [compare with Eq. (|T7|) ] , 

rr[{i>a}] = 

-WriUall+E / Ja[x)^a{x)-^S[^\{^,}]. 

(95) 

Its flow with the infrared scale k is described by an exact 
RG equation (ERGE))^iS which, after accounting for 
the curvature, reads 



atrf[{$a}] = ^ E / / (i-/3Mi)(i-/3M2) 

X dtUk-^a^a, (a , 3l2)<(i,.^)(a„.,) ii'^a}], 

(96) 

where t ~ log(A:/A); the (modified) propagator v\^'^ is 

defined as the inverse of F^ + TZk in the sense defined 
in Eq. (gHl), i.e., 

V / (1-2/3^3^3) Pf'^.^ w Jvf] w , 

+ 7?.fc; (03,^3) (02,^2)^ = + l30ldi)6aia2Sx^x^ , 

(97) 

(2) 

where F^. [{^"0}] is the second functional derivative of 
the effective average action with respect to the superfields 

and 4^3;^ = S^''-\xi - X2)5e^e^ [with, again, Jg^ Sg^g^ = 

(1 + pOiOi)]. Here and in the rest of this section, we 
omit the superscript (/3) in the functional derivatives of 
Ffe in order to avoid the awkward proliferation of super- 
scripts. (Note that conventional translational invariance 
being explicitly broken in the curved Grassmann space, 
Fourier transforming is of no use when (3^0.) 

Due to the properties of the infrared cutoff functions, 
the effective average action reduces to the standard effec- 
tive action F^'^Hi^a}] when k ^ 0. The initial condition 
at the microscopic (UV) scale, when — >■ A, is more 
subtle. By using the definition of the effective average 
action, Eq. ([M]). and of the regulator, Eqs. (|93l94p . and 
after a change of integration variables, one obtains the 
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expression: 



/n 
n2?XaCXp{-5(«[{$, + Xa}] 
a=l 

+ T. f ^klli'^- + Xa}]Xa{x) Rk{\x, - X2I 

/ Xa(a;i,0Xa(a;2,^) - ^ / -Rfc(|a;i - a;2|) 

n . n „ 

X(E Xa{x,A)){Y. Xa{x2,0))}- 
a=l Ji. a=l J I 

(98) 

If one requires that i?fe(|a;|) diverges for all x when fc A 
while stays bounded, the term in Ru in the above 

expression acts as a delta functional for all the superfield 
variables Xa- As a consequence, 



(99) 



i.e., the effective average action reduces to the bare action 
defined in Eq. 



B. Properties and role of the cutoff functions 



We first go back from Eq. (|92|) to the original formula- 
tion with the ifa fields in the presence of disorder (setting 
the fermionic sources to zero). The generating functional 
can then be expressed as 



Zf)[{j„jj](x 



det 



5ipa{x) 



Vh 



exp [ 



2(AB-i?,fc(g2))J 



det (As 

y)fa{y) - 



-Rk) 1=1 

h{x) - Ja{x) 



n y 



Rk exp / Ja{x)^a(x). 



(100) 

One can sec that i?^ is added to the second functional 
derivative Sq \ip\ of the bare action and plays the role of 
an infrared cutoff at scale k for the fluctuations of the Lp 
field, as in the NP-FRG of pure systemsi^ More specif- 
ically in the present case, a large enough cutoff function 
^ f 2) ^ 

i?fe ensures that the operator [(p] -I- Rk is definite pos- 
itive which guarantees that the stochastic field equation 
has a unique solution. This is certainly true at the UV 
scale A where Rk^A diverges. At the beginning of the 
flow, the regularized theory is therefore "ultralocal" in 
the Grassmann subspace. Since as discussed in section V 
the "ultralocal" cumulants Wp and Fp are independent 
of the curvature /3 when uniqueness of the solution is en- 
forced, the regularized theory at the start of the RG flow 
is also invariant under the superrotations. 

The other cutoff function Rk on the other hand (with 
the requirement that it is always positive) reduces the 



variance of the random magnetic source, i. e. reduces the 
fluctuations of the bare disorder. From this, one gets 
the additional constraint at the beginning of the flow, 
when the random fleld distribution is not yet renormal- 
izcd, that As > i?A(g^) for all q's. Note that as one 
follows the RG flow by reducing the IR scale k, Eq. pOOp 
is still valid but rather useless: it is indeed more con- 
venient to work with renormalized quantities, namely a 
renormalized disorder and a renormalized action or effec- 
tive action. When k decreases, so docs the cutoff function 
Rk , which presumably leads at some point to a breaking 
of the "ultralocal" property when /? is finite. One should 
however keep in mind that the relevant solutions are as- 
sociated with a renormalized random functional and no 
longer with the bare action, so that the deviation from 
"Grassmannian ultralocality" may be small and, even 
further, vanish along the flow (see below). 

Finally, the superfield formalism offers a way to con- 
strain the cutoff functions by relating them. As noted 
before, the action in Eq. (|42p is invariant under the su- 
perrotations when the supcrsources are taken as uniform 
in the Grassmann subspace (i.e., Jb = Kb = Kb = 0) 
for all copies but one and when the curvature /? is set to 
zero. The regulator ASk can be made explicitly invari- 
ant under the same conditions by choosing TZk,aa to be a 
fimction of the super-Laplacian Ass only. As a result)^ 



-ABd„2Rk{q^), 



(101) 



where q denotes the momentum in Euclidean d- 
dimensional space. Note that the above expression in- 
volves the strength of the bare disorder whereas, as also 
discussed above, a more general relation should allow for 
a proportionality factor expressed in terms of renormal- 
ized quantities rather than bare ones. This will be dis- 
cussed in section VII-B. 



C. ERGE for the cumulants 



Aiming at deriving exact ERG flow equations for the 
"cumulants" ^\^J , we first rewrite more explicitly the 
ERGE for the effective average action. 



a^rf = ! / (dMq'-) I (i-2/?0i0i)E 
^Si..-,,,)(a.,.^.)[i'f''}] + 5*i?.(g^)^ ^ {i~Pho,) 

(l-/3M2)E^SL-,.,)(a..,,,)[{<J'4]}, 
0102 ^ 

(102) 

and we use the expansion in increasing number of sums 
over copies. For the expansion of the full propagator 
P^*^' appearing in the right-hand side, it is convenient 
to introduce notations as follows. A generic matrix 
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Aai3Li)(o2i£2) ^'^ decomposed as 

AiaJ{$a}] = Sa.a^Aadi^a}] + Aa^aAi'^a}], (103) 

where we have dropped the exphcit dependence on the 
coordinates. In the above expression, it is understood 
that the second term Aaj^a2 no longer contains any ex- 
phcit Kronecker symboL Each component can now be 
expanded in increasing number of sums over copies, 

Aa, = ^["1 [$aj + ['^'-1 I*"] + • • • (104) 

6 

(105) 

where the superscripts in square brackets denote the or- 
der in the expansion (and should not be confused with 
superscripts in parentheses indicating functional deriva- 
tives). The A'^l's and ^[^I's are symmetric under any 
permutation of their arguments, and their functional 
form is independent of the number n of copies (taken 
as arbitrarily large). 

After inserting Eq. ((50)) for F^''-' and the above ex- 
pressions applied to the propagator 'pjf\ we obtain a 

hierarchy of ERGE's for the cumulants r^^-'. The first 
equations read: 

/ / (l-/30-i^?i)(l-/3^"2^2)Pl°l_,,^)(,,^)[<i>i]|, 

(106) 

dtRkiq") I I {I ~ i3W{i - (3e2e2)x 

-I- perTO(12) j-, 

(107) 

and so on, where perm{\2) denotes the expression ob- 
tained by permuting $i and 4>2; for clarity we have omit- 
ted the superscript (/?) in the right-hand side. 

The components of the propagator, and V]^ , can 
be expressed in terms of second derivatives of the cu- 
mulants Vkq by inserting into Eq. (^7)) the expressions 
P04ll05p with A equal to the second functional deriva- 
tive of Eq. (j50p . The algebraic manipulations are given in 



Appendix IB] One finds for instance that the propagator 
T'l^' can be symbolically expressed as 

vt\^i] = {v^'^^i]+RkUy\ (108) 
which means that 

^ (1 - 2ph0z) ^E^,, (j'-kL,., + Ru;.,.2U0_,o}j 
= Ue^ej'^''\xi ~ X2) , 

(109) 

where we have introduced Ug^g^ = (1 -I- f30i9i)Sg_^g^, and 
that pI^' is given by 

= / / (l-2/30303)(l-2W4) 

^il.3[*i](rSi3,.J<i'i:*2] - {i + p~eMi + phe^) 

(110) 

Eqs. (|TU5|) and (|TTU)) can be inserted in Eq. pilS)) . which 

provides an ERGE for f^'^-' only expressed in terms of 
cumulants associated with the effective average action. 

After introducing the short-hand notation dt to indi- 
cate a derivative acting only on the cutoff functions (i.e., 

dt = dtRk S/SRk + dtRk S/SRk), Eq. pII7|) can now be 
rewritten as 

dA^2 [*1'*2] = m / / (1 - "^I^Wd - 2/?04e4)x 

^1°^, (rif^.. 'i>2] - ri:3;2.^,,.. ['J'l' *2]) 

+ ^ii^ i'^i^'^i] rgu,. [I'l, 2] + ipi:;^ $2] X 

(rgL,.3 *2] - (1 + poMi + PhoA;x,x,) 

+ perm{\2) j^, 

(111) 

where, again, perm(l2) denotes the expression obtained 
by permuting $i and <f>2- Similar ERGE's are obtained 
for the higher-order cumulants. A graphical representa- 
tion of the hierachy of ERGE's is provided is AppendixlU] 
This provides a hierachy of exact RG equations for 
the cumulants of the renormalized disorder (including 
the first one which leads to a description of the ther- 
modynamics). One should note that (i) the cumulants 
arc functionals of the superfields and contain full infor- 
mation on the complete set of IPI correlation functions, 
(ii) the flow equations are coupled, the (p+l)th cumulant 
appearing in the right-hand side of the equation for the 
pth cumulant, and (iii) to obtain the flow equation for 
rfcp[$i, $p] with its full functional dependence on the 
p field arguments, one needs to consider at least p copies 
of the original system. Formally, the whole hierarchy of 
flow equations for the cumulants can thus be obtained by 
considering an arbitrary large number of copies. 



16 



D. ERGE for the cumulants under the hypothesis 
of "Grassmannian ultralocality" 

As such the hierarchy of ERGE's obtained above is 
expressed in terms of superfields and superspace coordi- 
nates and is untractable in general It is instructive to 
consider the simphfication that arises under the hypothe- 
sis of "Grassmannian ultralocahty" . This is achieved by 
considering the flow equations when the cumulants are 
evaluated for "physical" configurations $a(ai) = (t>a{x) 
of the superfields, i.e. configurations that are uni- 
form in the Grassmann subspace. Then, from Eq. ([7n|) . 

rg^ .., $p] = rg^ [01, Wc next insert the 
results of section V in the ERGE's and we introduce the 
following "hat" and "tilde" propagators, 

pi%] = {^^S[<p]+Rky' (112) 

and 

pfi[0i,<^2] = pf'[</'i](ri^/^[</'i,02]-i?fe)pfi[02], (113) 

which are obtained from Eqs. (|108mi0|) with, as in 
Eqs. = (l + PWSe^e.Pi'hcj)] and 

^k^ej<l>uh] = (l + /3^iei)(l + /3^2^2)Pf'[0i>2]. 

The exact flow equation for the first cumulant of the 
rcnormalizcd disorder is finally obtained as 

dtTki = 

(114) 

where the explicit dependence on the Euclidean coordi- 
nates has been momentarily reinstalled. One similarly 
derives an ERGE for the second cumulant, 

df^k2 [01,</>2] = 

JX3X4 

X (rgl,.. [<^i>'^2] - Rk;x,x,) +perm{12)}, 

(115) 

where perm{12) denotes the expression obtained by per- 
muting and (j)2. 

Generically, the flow of F^p (/)p] involves three 

types of quantities: the propagators P^°' and pj^\ second 
functional derivatives of T^p in which all the arguments 
are different, and second functional derivatives of T^p^i-^ 
with two of their arguments equal to each other (for a 
graphical representation of the hierachy of ERGE's, see 
Appendix [C|) . Wc will come back in more detail to the 
structure of these flow equations in the following paper, 
but we note for future reference that these equations are 



independent of the auxialiary parameter [3, and so is their 
solution if the initial condition is itself "ultralocal" . 

Finally, wc point out that the above ERGE's coin- 
cide with those previously derived without the supcrficld 
formalism by means of an expansion in number of free 
replica sums, when evaluated at T = O j^^-^^ The same 
is true for the ERGE for all higher-order cumulants. It 
is however important to stress that our previous replica 
approach provides no obvious way to make the superro- 
tational invariance (or lack of it) explicit and, therefore, 
neither to relate the two infrared cutoff functions and 
Rk nor to provide guidance for devising nonperturba- 
tive approximations to the ERGE's that do not explicitly 
break the underlying supersymmetry (see the companion 
paper) . 



VII. GROUND-STATE DOMINANCE IN THE 
NP-FRG 

A. Taking the limit of infinite curvature (zero 
auxiliary temperature) 

We now come to a central step of our approach. Intro- 
ducing the auxiliary temperature /3~^ has allowed us to 
place the supcrficld formalism on a firm ground. How- 
ever, being ultimately interested in studying the ground- 
state properties of the system, we would like to take 
the limit /3^^ — > in the exact flow equations derived 
above. To make the dependence on /? explicit in the 
ERGE's, we flrst rescale the Grassmann coordinates and, 
accordingly, the auxiliary flelds: (8,0) ~ /3^^/-^(tZ;, w), 

(1/;,^) = /3i/2(^^^)^ and (j> = so that $(^) = $(w) 

with $(aj) ~ (j) + Cjij] + V'^ + Cjijj4>. 

The cumulants associated with the effective average 
action can be formally written as 




(116) 



where in the right-hand side, as in Eqs. (|44I45I) . q E 
{1, denotes 

{'J'a,(^,),(l + fM«)5e,<f>a,(^,),A,^$,^(0,)}. (117) 

When changing Grassmann coordinates from 9 to w, we 
assume that Fg^[l, ...,p] = Fg^[i, where 1 denotes 

{^ai(wi), (1 + 5WiWi)a^^$ai(wi), Aj^^$ai(tiii)}, and so 
on, with the same functional form for Fg'' . This guaran- 
tees that the derivatives of the "non-ultralocal" contri- 
butions to the F^-* 's do not come with increasing factors 
of (3 which would completely spoil the limit of inflnite 
/?. We shall support this argument by studying the first 
"non-ultralocal" corrections (see below). 
Then, taking into account that an integral over a 
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Grassmann variable acts like a derivative leads to 

ri?[*a„...,$aJ = /3^ / ... / rg^[i,...,p]. (118) 



In addition, by using the identity S/S^ai^) = S/S^ail^), 
which can be proven by considering the components of 
the superfield in the old and the new coordinate system, 
we arrive at 
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<[l,-.,;5], 



(119) 

where the functional dependence on the Euclidean coor- 
dinates is left implicit. 

After the above preliminaries, we consider the ERGE's 
for the cumulants, e.g. Eqs. (|106lllip for the first two cu- 
mulants. We have to study the explicit /3 dependence of 
the propagators 7^1°' and P^"' . After changing the Grass- 
mann coordinates in Eqs. (|109lll0p and using Eq. pi9p . 
we easily obtain 



(120) 



where Tr indicates a trace over the Euclidean momenta 
(which are not shown explicitly) . Similar expressions are 
obtained for the higher-order cumulants. One actually 
finds a structure that is analogous to that derived in our 
previous replica approachpi^ when considering the (bath) 
temperature T: the cumulant of order p comes with a 
factor T~P, its qth derivative with a factor T~'-p~'^\ the 
"hat" propagator has a factor 1/T and the "tilde" prop- 
agator no explicit factor of T. The limit /3 — > od allows 
one to drop the terms in the above flow equations that 
have an explicit factor of 1//3, namely the first term of 
the right-hand sides. It should however be kept in mind 
that as in the "thermal" case when T — > 0, the limit 
1//3 — > is expected to be nonuniform in the (super)field 
dependence. This will be discussed further down and in 
the following paper. 

Having introduced the change of Grassmann coordi- 
nates and the formal way to study the limit /3 — > oo, we 
can go one step beyond in the analysis of the ERGE's. 
To do so, we consider superfields that are uniform in the 
Grassmann subspace, i.e. $a = = (j^a, and we gen- 
eralize the results of sections V-B and VTB to the case 
where the cumulants vfj have a generic "non-ultralocal" 
component. 



rij^^[<i>i(aii),...,$p(c^,)] 



and 



*i>'f2]-^l:i^^J<i'i,<i'2], (121) 



where <& is the superfield defined with the new coordi- 
nate system and the rescaling of the auxiliary fields (see 
above) . 

The ERGE's for the cumulants can then be reexpresscd 

as 

' ^11,., + ^i;L,., ^-i]) + dtRky- (122) 



(l-cSlC.l)(l-<i2a.2)^il,„J<i>l]l>, 



d,l I C[i,2]=ia,Tr 



(1 - 2W3W3)X 



(1 — 2w4a;4) 



-(20) 



$1 



r(llO) 



$1 



-(20) 



$1,$2 



$1,$2 



-V 



[0] 



$1, $2 



-(11) 



k2:uj_. ,LJ 



$i,$i 

$1,$2 



{1 + uj4LUi){l + uj3UJ3)Rk] +perm{12) 



(123) 



-,{P)NUL 
kp 



[l,...,p\ 



(124) 

involves derivatives of the superfields in 
the Grassmann directions and is equal to zero when the 
superfields are uniform in the Grassmann subspace: as a 



where T 



(I3)NUL 
kp 



,pp(/})C/Lr 



. The second 



result, r^^^[^i, (j)p\ = pPT^^ wi,--, vp\ 
functional derivative of the effective average action can 
be decomposed as in Eq. (|103p and, by applying the WT 
identities associated with invariance under the isometrics 
of the curved Grassmann subspace (see section IV-D and 
Appendix |D]) , one finds the following general structure 
for the "hat" and "tilde" components: 



(1 + UJlUJl)5^^ui^Ak-ai [{4>a}] 
-f (1-f UJ1UJ2 + Ul2(^l)Bk;ai [{4>a}] 



(125) 



and 

(126) 

where we recall that a factor 1//3 is present when chang- 

^(2) 

ing variables from to w in . 

In addition, it can be proven that Ak-ai and Ck-aia2 
are obtained only from the "ultralocal" part of the cu- 
mulants, namely (leaving again implicit the dependence 
on the Euclidean coordinates). 



[</'ai , (pa 



(127) 
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Ck;aia2[{'Pa}] 



^UL{11) 
k2 



C/L(110) 
k3 



(128) 

whereas B^ an whose precise expression is not partic- 
ularly illuminating at this point, involves the "non- 
ultralocal" part of the cumulants. As a result, Bk-ai is 
equal to zero when the effective average action is purely 
"ultralocal" in the Grassmann subspace. The demon- 
stration of the above property and of Eqs. (|125m28|) is 
provided in Appendix [Dl 

The full propagator 'Pk-(aiui_^)(a2ui^)-, which is the inverse 

of r^^^ + TZk, has the same structure as in Eqs. (|125ll26p 
with 



[{4>a}] = (1 + ^l^^l)Sui iWa}] 

+ (1 + OJ1UJ2 + '^2Wi)Qfc;ai [{0q}] 



(129) 



and 



/ A;;(aicj-|^)(a2^2) [{4>a}] = {l+U!lUJi)(l+UJ2^2)Pk;aia.,[{4>a}], 

(130) 

where Vk-a a-nd Vk;aia2 ^re expressed only in terms of 
"ultralocal" cumulants whereas Qk;ai also involves "non- 
ultralocal" terms and therefore vanishes when the latter 
go to zero. All these quantities can be expanded in in- 
creasing number of sums over copies. The zeroth-order 

components of the propagator, T'j"', Pj^^ and Q^^\ are 
then expressed as 



[^]={T]!I'^'\^]+Rk 



(131) 



which denotes an inversion in the sense of operators in 
Euclidean space. 



P. 



[0]r 



[Ojr 



i](r 



(132) 



and 

(133) 

where the Euclidean indices are omitted for simplicity. 
Note the (expected) correspondence between the above 

expressions for pj^^ and P|°' and those derived under the 
hypothesis of "Grassmannian ultralocality" , Eqs. (fTT2|) 

and dnni). 

Reinstalling the explicit dependence on the Euclidean 
momenta, we find the following exact RG flow equations 
for the "ultralocal" parts of the cumulants. 



p[0] r , 1 / T.C/L(11) 
^k:qq'l'Pl\\^ k2:q.q' 

^ -'qq' 



5^'^)[q-c^)Rk{q^))+^^ J dMci')Qtlqq[<l>ll 



qq' 



(rS"' - {2^f5^'\q - q')Rk{q^)) + perm{l2) 

+ ^^*/ , {^S^' ["^il^™^ ['^i''^^] +per™(12)|, 

(135) 

and similarly for the higher orders. 

Setting 1//3 = in the above ERGE's allows one to get 
rid of the terms that involve the "non-ultralocal" contri- 
butions to the efi'ective average action. The /3 — > 00 
limit therefore coincides with the RG equations for the 
cumulants obtained under the hypothesis of "Grassman- 
nian ultralocality" (compare with section VI-D). This 
shows that the ERGE's derived under the assumption of 
"Grassmannian ultralocality" describe the renormaliza- 
tion of the RFIM at equilibrium with a proper selection 
of the ground state. This is an important piece in the res- 
olution of the problem associated with the long-distance 
physics of the model. 



B. Illustration of the RG flow for "non-ultralocal" 
contributions 

We now illustrate the structure of the ERG flow for the 
"non-ultralocal" components of the cumulants by looking 
at the lowest-order correction to the first cumulant, as- 
suming that all other cumulants are purely "ultralocal" . 
This allows us to verify that this flow is well-behaved, 
thereby giving a direct confirmation to the arguments 
used in the previous subsection. More specifically, we 
consider 

Vkim = P (r^kilM^,)] 

-1 +in[$i(^i)](l+cSic^i)9^,$i(aii)9^,$i(aii)Y 



and for p > 2, 



(136) 



r,.,[$i,...,$p] = /3f / ... / r^/^[$i(^o,...,ci>p(wp)], 

(137) 

where we have omitted the superscript (/3). 

The ERG equations for the "ultralocal" components 
are the same as in Eqs. (|134ll35p and the propagators 

-Pf \<t>l Pk^ [4>i,4>2] and qf [(/)] are given in eqs. (il3U1133|) 
with B^"' [0] now simply equal to 



Bfm = -Ykm. 



(138) 



(134) 



The most convenient way to obtain the flow of Yk[(t)\ is 
to differentiate twice the ERGE for f^i in Eq. (|122p and 
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evalutate the resulting expression for a superfield config- 
uration that is uniform in the Grassmann subspace (i.e. 
$ = (/)). The final RG equation for Yfc [(/)] reads 



dtYk [0] = atTrjpf 1 



[Qr [</>] - ^1"' [^] ] ) + 2 [Q^ m - [^] ] ( pi:' ^] X 



1 



(2) 



(1) 



[Olr 



5(0] 



5(0] 



kl 



[0] 



fc2 



(1) 



5(0] 



^UL{3} 
kl 



^UL{22) 
k2 



2Pf 1 [0] X 



UL(21) 
k2 



2 



(2) 



(1) 



[4 



[0] 



k 

m]^ki 



ULCi) 



2 j,UL[3) 
^ kl 



(139) 

where Tr denotes a trace over Euclidean momenta (which 
are left implicit in the right-hand side of the equation) . 

Provided that Yk[(t)\ converges to a finite (nondiverg- 
ing) fixed-point value, its contribution to the fiow of the 
"ultralocal" components of the cumulants, which appears 
through the propagator in Eqs. (|134p and (jlSSp . 

can indeed be neglected when /3 — > oo. This will be 
shown in the companion paper. (More generally, it can 
be shown that the contribution becomes subdominant at 
long distance, i.e. when fc — >■ 0, even when j3 is large but 
finite.) 



C. Scaling dimensions near a zero-temperature 
fixed point and asymptotic dominance of the ground 

state 

To search for the fixed point that controls the criti- 
cal behavior of the RFIM associated with the sponta- 
neous breaking of the (global) Z2 symmetry, the NP- 
FRG fiow equations must first be recast in a scaled form. 
This can be done by introducing appropriate scaling 
dimensions J^ii^ Near a zero-temperature fixed point, it 
is convenient to introduce a "renormalized temperature" . 
As shown in Refs. [l6lll8| . this can be done by consider- 
ing (i) the strength of the renormalized random field , 
which can be defined from the vertex T^^2^ evaluated for 
a specific configuration of the fields and reduces to the 
bare value As at the UV scale A, and (ii) the amplitude 

of the field renormalization constant Zk, which is as usual 
(2) 

obtained from F^-^ for a specific field configuration and 
is equal to 1 at the UV scale. Specifically, 



Zk^ 



(140) 



An associated running exponent 0^ is defined from 

ek=dt\ogTk (141) 

whereas the (running) anomalous dimension 77^ is ob- 
tained as 

rjk = -dt\ogZk. (142) 

One may also introduce a running exponent 77^ = 2 — 
dk + m from 

-2r/fe =9tlogAfe. (143) 

Note that the two anomalous dimensions "q and 77 describe 
the spatial decay of the pair "connected" and "discon- 
nected" correlations functions (see section [V]) at critical- 
ity; this translates into 



Pl\q-^) 



-(4-))) 



(144) 



with r] <fj < 2r^i^. 

A systematic way to proceed is to introduce, on top of 
the canonical scaling dimensions for the effective average 
action and the related cumulants and of a rescaling of the 
Euclidean momenta and coordinates (q = q/k, x = kx), 
a rescaling of the Grassmann coordinates via a renormal- 
ized curvature (3k oc 1/Tk, with Tk defined in Eq. (|140p : 



where 



Pk = P 



Zk{k/Ay 



(145) 



(146) 



The renormalized curvature reduces to the (bare) curva- 
ture /3 at the UV scale A and diverges as fc"^ in the IR. 
It is worth pointing out that a symmetry between the 
rescaling of the Euclidean and of the Grassmann coordi- 
nates exists in the case where the dimensional reduction 
applies: then, = 2 so that all coordinates are rescaled 
by the same factor k ~ k^^^. 

We can now introduce dimensionless quantities. The 
dimensionless superfield is defined through 



<^{x,e) 



Pkk' 



d~2 



PZk 



1/2 



$ren(i,£), (147) 



which implies for the components: 



4>{x) 



PZk 



■„d-2 \ 1/2 



ZkPk 



(fix), 



(148) 



(149) 



etc. 



20 



The "ultralocal" component of the pih cumulant is 
rescaled as 

r^/^[0i,...,^J = fc'^(/3fcf 7.';'^[^i,-,¥'d, (150) 

and the "non-uhralocal" cumulant can similarly be put 
in a dimcnsionless form . For instance, the term B^^"^ [(f>] 

appearing in the definition of the propagator Q^^^ [ip] is 
expressed as 

B^°\(l>] ^ PkZkk^y]. (151) 

The dimcnsionless quantities will be systematically de- 
noted by lower-case letters (except for the supcrficld and 
for the coordinates, for obvious reasons). 

In addition, and following the discussion in section VI- 
B, the cutoff functions arc chosen according to 

Rk{q^) ^ Zkk^sif), (152) 

with s{x) such that s(0) > and s{x oo) = (see the 
companion paper and Refs. [l6lll8j ). and 

Rkiq') = 'Aks'if), (153) 

with s'{x) the derivative of s{x); the above form of 
the regulator then satisfies the superrotational invariance 
whenever A^/Zk ~ Ag. 

We focus here on the "ultralocal" components of the 
cumulants. The ERGE's for the latter can then be put in 
a scaled form which has the following generic structure: 

Qtlkv ['/'I ' ■ • • ' V'p] = -P(2 + r/fc - T]k)-ikp [V'l , ■ . • , V'p] + 

Pfc 

(154) 

where is a dimcnsionless beta-functional expressed 
only in terms of (dimensionless) "ultralocal" components 
and T^^^ is a dimensionless beta-functional that con- 
tains (dimensionless) "non-ultralocal" components. Ex- 
plicit expressions for the first cumulants are easily derived 
from Eqs. (|134ll35p and will be given in the following pa- 
per. Here, only the general structure is needed. 

Fixed points are then found by setting the left-hand 
side of Eq. (|154p to zero. This will be studied in detail in 
the companion paper. At this point we would only like 
to stress that the "non-ultralocal" contributions to the 
flow of the dimensionless renormalizcd cumulants come 
with a factor l//3fe. Even if /3^^ is not taken equal to 
zero, and provided the "non-ultralocal" contributions re- 
main bounded (see following paper), the flow leads to 
a fixed point characterized by the property of "Grass- 
mannian ultralocality" : one indeed expects l//3fe oc 
to flow to zero, i.e. Ok^o > 0, as already shown in com- 
puter studies^ and in Refs. [iTlflQj . Ground-state dom- 
inance is thus found asymptotically as the flow goes to 



the zero-temperature fixed point. Guided by our pre- 
vious workiiiiS and by that of Balents, Ledoussal and 
coworkers on the random manifold model;^'^'^'^ we 
anticipate that when the fixed point is characterized by 
a nonanalytic, cusp-like, dependence of the cumulants of 
the renormalizcd random field, the approach to the /c — > 
limit (which coincides to the fixed point obtained by set- 
ting rightaway /3~^ = in the ERGE's) involves a bound- 
ary layer, generically in \ipa ~ fbl/Tk^ This boundary 
layer is physically associated with the presence of rare, 
power-law distributed, "droplet" excitations above the 
ground state at the running scale fci^ 



VIII. CONCLUSION 

In this paper we have extended our nonperturbative 
ERG approach of disordered systems, which was pre- 
sented in the previous articles of this series J^ii^ The ob- 
jective was to discuss the property of dimensional reduc- 
tion and its breakdown in the RFIM from the Parisi- 
Sourlas^ perspective of an underlying supersymmetry 
and its breaking. To this end, we have reformulated the 
ERG in a supcrficld formalism. 

We have identified two sources of problems in the 
Parisi-Sourlas supersymmctric formalism. One stems 
from the presence of metastable states which, due to the 
resulting multiplicity of solutions to the stochastic field 
equation that describes the long-distance physics of the 
RFIM, prevents the selection of the ground state in the 
random generating functional. The other one was already 
discussed in our previous papers and comes from the fact 
that a single copy of the original system is considered. As 
such, the formalism is therefore unable to describe rare 
events, such as avalanches and droplets, that manifest 
themselves as nonanalyticities in the field dependence of 
the cumulants of the renormalizcd disorder. 

We have provided ways to cure both problems, through 
the introduction of a weighting factor with an auxiliary 
temperature and through the use of multiple copies. The 
resulting theory involves superfields in a curved super- 
space whose curvature is related to the inverse of the 
auxiliary temperature. The presence of metastable states 
leads to a breakdown of a formal property which we 
have called "Grassmannian ultralocality" and is asso- 
ciated with the fact that a unique solution is incorpo- 
rated in the random generating functional. On the other 
hand, as will be discussed in detail in the companion 
paper, nonanalyticities originating from the presence of 
avalanches in the ground state as one varies the applied 
source trigger a spontaneous breaking of supersymmetry, 
more precisely of "superrotational invariance" . 

Through the introduction of an appropriate infrared 
regulator which guarantees that the stochastic field equa- 
tion has a unique solution at the initial condition of the 
fiow and which can be chosen to satisfy all (super) sym- 
metries of the theory, we have derived exact RG flow 
equations for the cumulants of the renormalizcd disor- 
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der in the effective average action formalism. We have 
shown that in the Umit of infinite curvature, i.e. of zero 
auxihary temperature, the hierarchy of exact FRG equa- 
tions coincides with that obtained under the hypothesis 
of "Grassmannian ultralocality" . Through this proce- 
dure, the ground state is properly selected and is the 
only solution that contributes to the random generating 
functional. We have moreover found that the corrections 
to "ultralocality" that are present for a finite curvature 
become subdominant at long distance as one approaches 
the expected zero-temperature fixed point. 

In the following paper, we investigate the exact FRG 
equations derived for the cumulants in the limit of infi- 
nite curvature and show that nonanalytic behavior of the 
effective average action in its field arguments is related to 
the appearance of a spontaneous breaking of the super- 
rotational invariance. We also devise and study a non- 
perturbative approximation scheme that does not break 
explicitly the superrotational invariance and allows us to 
study the critical behavior of the RFIM and the break- 
down of dimensional reduction as a function of the spatial 
dimension d. 



Appendix A: "Metastable" states and 
nonanalyticities in the zero-dimensional RFIM 

In section lllli we have introduced the property of 
"Grassmannian ultralocality" by considering the simple 
case of the d = RFIM at zero temperature. We use 
this same toy model here to illustrate on the one hand 
the different proposed ways to construct generating func- 
tionals when "metastable" states are present and on the 
other hand the nonanalyticities that may be present in 
the field (or source) dependence of these functionals and 
of the associated Green's functions. In the d = RFIM, 
a variable ip is submitted to an interaction characterized 
by the action: 

S{^)^UB{ip)-{J + h)ip, (Al) 
with, as in Eq. ([T]), 

Ub{^)^1^' + ^_^' (A2) 

and the random field is Gaussian distributed with vari- 
ance Ab- 

In the case where r is negative, the "stochastic field 
equation" which reduces to Ug{(p) = J + h has three 
branches of solutions: 

(1) (p-iJ + h), which exists for J + h < Jc with Jc = 
2M /ijTF 

3 V " ' 

(2) (/)+(J + h), which exists for J + h > — Jc, and 

(3) (j)f{J + h), which exists for — Jc < J + h < Jc. 

In the region —JcKJ + hKJc, the equation has then 
three solutions, two minima, 0_ and (^_|_, and a maxi- 
mum, (j)^(J+h). This is a situation where one encounters 
multiple solutions or "metastable states" . 



We focus on a "Green's function" which, when a sin- 
gle solution is present and all constructions lead to the 
same results, is obtained from the second derivative of the 
second cumulant of generating function W(J, J), namely 
[compare with Eq. ([7])], 

Mjr+W¥{jT+h)) 

-{ip{Ji+h)) {ip{J2 + h)). 

At T = 0, the thermal average (ip) reduces to the ground 
state (j)o, i.e. the absolute minimum of the action. In 
the Parisi-Sourlas supersymmetric construction, (ip) is 
replaced by </)+ -I- 0_ — (/)| (in the range where the three 
solutions coexit) and in the construction proposed in this 
work (see section III-B), (ip) is replaced by a Boltzmann- 
like weighted average with an auxiliary temperature 

As we want to study the possible nonanalyticities in 
the dependence on the sources Ji , J2 as J2 — >■ Ji , we 
rewrite Ji = J — SJ, J2 = J + SJ and consider the limit 
(5J — > 0. Moreover, we consider the difference between 
the Green's functions 

MJ + SJ + h)(l)^{J - SJ + h) - JjJTW, (A4) 

where we have dropped the difference in the disconnected 
pieces that is easily shown to be at least of 0{SJ^) and 
0* is given by any of the three above definitions (ground 
state, Parisi-Sourlas weighting, Boltzmann-like weight- 
ing). The expression in Eq. (jA4| is symmetric in the 
inversion SJ — >■ —SJ and so, without loss of generality, 
we choose SJ>Om the following. 

a. Ground state 

If 0* is taken as the gound state, it then presents a 
discontinuity at J + h = (being equal to (f>- when J + 
h < and to 0+ when J + h > Q). It is easy to show that 
0*(J + SJ + h)(j)^{J - SJ + h) - <f)^{J + hY is of order 
0{SJ^) everywhere except in a region centered around 
J + h = Q and of width 2(5 J, where the function is of 
order 1 (actually — 12|r|/u). Averaging over the disorder 
therefore generates a linear cusp in SJ , 

(f)^{J + SJ + h)<i)^{J - SJ + h) - <i)^{J + hY = 

24 P (A5) 

- exp(--— )|M + 0{SJ'). 

v27ritAs ^i^B 

h. Parisi-Sourlas weighting 

If we use the prescription that 0* = -f 0- — (/"i when 
the three solutions coexist {i.e. for — Jc < J+h < Jc) and 
0* is given by the unique solution otherwise, the function 
4>*{J + h) is now continuous but presents a nonanalytic 
behavior around ± Jc. When approaching Jc from below, 
one has (/)*(J) ~ (l)+{Jc) - A^yjc - J + 0{J ~ Jc) with A 
a positive expression that is not worth to be given here. 
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A similar behavior is observed when J approaches —Jc 
from above. It is now easy to show that (/)++(/)_— goes 
to zero as 0(5J^) everywhere except in regions of width 
5 J around J + h = ±Jc. In these regions, the function 
behaves as 0{\/6J). After averaging over disorder, one 
therefore obtains a singular behavior 



MJ + SJ + h)(l),{J - SJ + h) - (l)4J + hp cx |<5Jp/2^ 

(A6) 

as 5 J — > 0. However, this nonanalytic behavior originates 
from branches that arc not absolute minima of the action 
and should therefore be discarded in studying equilibrium 
behavior. It is an artifact of the Parisi-Sourlas construc- 
tion. 



c. Boltzmann-like weighting 



In this case, one has 



)/(e- 



-/3S+ 



-/3S_ 



e-'^'^t) where 5*+ = S{(j)+), 



etc, when the three solutions coexist and 0* is given by 
the unique solution otherwise. As in the Parisi-Sourlas 
weighting, (/)*(J -I- h) is continuous with a nonanalytic 
behavior around ztJ^. In this case however, the result- 
ing nonanalyticity in |(5Jp/^ in the difference of Green's 
functions, Eq. (jA4p . is exponentially damped by a factor 
exp(— /JAS') where A5' is the difference of action between 
the ground state and the first excited state. 

The behavior of around J -\-h = Q is continuous and 
regular, but with a rapid change of 0(1) over a narrow 
interval of width 0(/3"^). After averaging over disorder, 
this gives rise to a "boundary layer" in /3|(5J|, i.e.. 



0* {J + 5J + h)(j)^ {J - SJ + h) - 0* (J + 



(A7) 



with f{J,y) = /2(J)y^ H when y and /(J,y) 

fc{J)y when 2/ — > oo. 

When (3 — > oo, the spurious nonanalyticity in |(SJ|'^/^ 
vanishes exponentially fast whereas the boundary layer 
leads to a linear cusp in |<5J|. One exactly recovers 
ground-state dominance [and it is easy to check that 



fciJ) - 
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{V2ttuAb 



■exp(-JV(2AB))] 



Appendix B: Algebraic manipulations for the 
expansion in number of sums over copies 

In this appendix wc provide a brief survey of some of 
the algebraic manipulations that are associated with the 
expansion in increasing number of sums over copies. We 
recall that the pth term of such an expansion involves 
exactly p copies and is a symmetric and continuous func- 
tional of either p copy superfields or p copy supersources 
(depending on whether one works with superfields or su- 
persources as fundamental variables) , the functional form 
being independent of the total number of copies (which 



however must be larger than p to give access to the full 
functional dependence) . 

For instance, to relate the terms of the expansion of 
the effective action r[{$a}] [wc omit the superscript (/3) 
throughout this appendix] to those of the expansion of 
M^[{j7a}], i.e. the cumulant of Wh, one can use a pedes- 
trian method that starts with the Legendre transform in 
Eq. (|Tf)) . Keeping the superfields {$o} as variables, we 
also expand -yV[{ Ja}] + J2a L '^a{x)^a{x) in number 
of sums over copies as 



$^(-Wi[X[{$;}]]+^J,(0;{$/})cI>,(^) 



n n 



(Bl) 



j3>2 ai—1 ap — 1 



where the supersourcc has the following expansion: 



{i + l38e)Ja{0;{'i>f}) = r[''l,[^a 



EE- 

P>1 CLl — l 



{-l)P 



E^(p + 1)! 



-(10. .0) 

p+i.e.. 



[$a,$ai,-,$aj 



-(0..01) 
p+l,..( 



[$ai,-,$ap,$a 



(B2) 

Here, we have dropped the explicit mention of the depen- 
dence on the Euclidean coordinates to keep only that on 
the Grassmann ones. 

We next proceed to a term-by-term identification of 
both sides of Eq. (|47| and we symmetrize the output if 
necessary. The first orders are easily derived and are 
given in Eqs. (|5m55p . Higher-order terms are similarly 
derived. One finds that the pth. term Vp[^ai, --j^a ] is 
equal to (-l)PWp[ J[$aJ, -^Ji^a,]]: where J[^a] is" the 
nonrandom supersourcc defined in Eq. (|54p . plus terms 
involving cumulants Wq with q < p that are also func- 
tionals of the nonrandom supersourcc. 

We also have to consider the expansion in increasing 
number of sums over copies of matrices, as in section 
VI-C. A generic mat rix ^(aio; j(a^ j;^) [{$a}] can then be 
expressed as in Eqs. (|103ll04ll05p . Algebraic manipula- 
tions on such matrices can again be performed through 
a term-by-term identification of the orders of the expan- 
sions. For instance, the inverse B = of the matrix 
A in the sense of Eq. (|49)) . i.e. 



(1 -f pe,e,)Sa,a, 5^,x, = E / (1 - 2/3^3^3) 



as -3 



(B3) 



X B, 



(ai^j)(a3^3) 



[{1>a}], 



can also be put in the form of Eq. (|103p and its compo- 
nents, Bai and Baia2 1 can be expanded in number of sums 
over copies. The term- by-term identification of both sides 
of Eq. (jBSp leads to a unique expression of the various 
orders, ^Wand B^\ of the expansion of B in terms of 
the ^['I's and ^^'I's with q < p. The algebra becomes 
rapidly tedious, but the first few terms are easily derived. 
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^[01 is the inverse of ^["1 , i.e., 

(1 - 2/30"303)^|'1.3[$i]-4:|'k[*iJ = (1 + 

(B4) 

whereas 

^^2] = - / / (1 - 2PM3){1 ~ 2/30404) 



and 

sW[$i|$2]=-/ / (l-2/30303)(l- 2/30404) 



(B5) 



(B6) 



6W[$i,$2|*3] = - / / (l- 2/30303)(l- 2/30404) 

+ ^W,J$l|$3]^|'l,J$i,$2] 

(B7) 

etc. These equations allow us to relate the expansion of 
the full propagator to that of the effective average action 
and to derive the ERGE's for the cumulants in section 
VI-C. 



Appendix C: A graphical representation of the 
hierarchy of ERGE's for the cumulants 

In this appendix we give a graphical representation of 
the flow equations for the cumulants T^p. We concen- 
trate on the replica structure and remain general on the 
way the momentum and Grassmann coordinates must 
be integrated over in the loops. Similarly, we use a nota- 
tion with "simple" fields but superfields could be used as 
well. In doing so, we obtain results which apply to the 
curved space, to the "ultralocal" case (when curvature 
drops out), but also in the replica formalism of papers I 
and 11^^ 

It is convenient to represent the first cumulant by a 
filled circle 



1 • 



rfei[</'i], 



(Cl) 



the second cumulant by two filled circles joined by a 
dashed line 



2 • 

1 * 



fe2 <Pl,<P2 



(C2) 



the third cumulant by three filled circles joined by dashed 
lines 



,«3 

Y 

1 • 



(C3) 



etc. The rule is that each circle corresponds to one copy 
(or replica), i.e., to one field 0a. We represent the zcroth- 
ordcr "hat" propagator by a line: 

^ K^^M- (C4) 

The vertices are represented by extracting legs from 
the dries. For the first cumulant, the graphical rule is as 
usual 




±3, 



kl.x-. ■■■x^ 



(C5) 



For the second cumulant, we must specify on which 
circle we extract the leg. We represent the vertex 
^'k2'^--x^ .y -v with qi legs in the circle asso- 

ciated with the copy field 0i and q2 legs in the circle 
associated with the copy field (1)2 ■ Graphically: 

^1 



y.1 



^(^91,92) 



y 



■12 



(C6) 

Some caution must be taken in the particular case of 
r^2'^^ as it always comes with the cutoff function Rk. 
We therefore use a special graphical rule in this case: 

^ rl^2i,.J'^i:^2] 

- (1 + /30i0i)(l + /30202)i?fc(a; - y). 

(C7) 

Finally, when needed, we encode the time derivative of 
the regulating functions in the following way: 



X 



dtR 



X 



X 1 



-dtR 



(C8) 



(C9) 



Let us now derive the graphical representation for the 
inverse of F^^) + TZ. To this end, we first write F^^) + TZ 
as 



(CIO) 
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In this expression, the external legs are amputated, 
i.e., do not include propagators. We can now use the 
inversion formula {A - = A^^ + A^^BA^^ + 

A-^BA-^BA-^ + ... graphically: 



and 



n 



Po 



r 



(Cll) 

In these diagrams the external legs now come with prop- 
agators. 

We now easily obtain the flow equation for F by con- 
tracting the previous expression with dtlZ. We readily 
find 




We are now in a position to write the flow equations for 
the cumulants by retaining the diagram with 1 copy, 2 
copies, etc, and appropriately symmetrizing the results. 
These expressions simplify if we rewrite them in terms of 
the dt operator. Recall that the operation dt corresponds 
to a derivation of the t dependence of the cutoff functions 
only. When performing this operation, we must use the 
rules 



dt 
dt 

This leads to 

dtri[<i>] 
2^ 

Y _ 

1 • 



(C13) 



(CM) 



(C15) 



2 • 2 

2 



1^ 



^0 



dtT^ = ]^dt\ -2 



r 

6 



3 X -3»-^-« + 3 



(C17) 



(C16) 



where it is understood that the right-hand side of the 
above expression should be symmetrized with respect to 
the copy fields. Similar expressions can be derived for 
the higher cumulants. 

In the case where "ultralocality" is satisfied (or im- 
posed), the last diagram of Eq. (jClSp . the last two dia- 
grams of Eq. (|C16p , and the last two of Eq. (jC17p give 
no contribution because the propagator, which is propor- 
tional to Sg_^g^, is contracted with a single copy, i.e., with 
9 1 = 9_2, or because two propagators connect the same 
copy, giving a contribution ((5e^ei^)^ ~ 0. Note however 
that these diagrams give nonzero contributions if "non- 
ultralocal" terms are present in the action. 



Appendix D: "Non-ultralocal" contributions; 
Grassmann structure of the 2-point vertex functions 

We derive in this appendix the dependence of the 
2-point vertex functions on the Grassmann variables. 
We show in particular that the "ultralocal" and "non- 
ultralocal" parts of the effective average action lead to 
distinct Grassmannian structures. 



1. Consequences of the WT identities 

Let us start with r^^"*. We derive Eq. (|6T|) with re- 
spect to ^a{x,9_i) and ^b{y,d_2) and evaluate the result 
in a field configuration uniform in the Grassmann vari- 
ables, i.e. $e = (^e,Ve G {!,..., n}. Focusing on the 
part proportional to 6ab and observing that the term with 
Grassmannian derivatives of the field vanishes, we get the 
identity 

de, ((1 - Phe,)f^il^,,^+de, ((1 - Phe2)f'-:i^,._) = o. 

(Dl) 

By using the symmetry of the 2-point vertex function 
under permutation of the indices and the fact that it has 
a vanishing number of ghosts, we can parametrize Fq as 

fi;(e,a;i),(e2X2)['t'?^«}] ^a;xtx._ + Ba-^.x^A^lOl + 6262) 

+ Ca [9,92 + 920,) + D 



(D2) 
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where Aa:xiX2, and Da-xxx2 are function- 



als of the 0e's. Pluging this expression in Eq. (|Dip leads 
to a set of constraints, 



Da - PBa, 



(D3) 
(D4) 



where we have dropped the expHcit dependence on the 
Euchdean coordinates. From the solution of the above 
equations, we derive the following form of the 2-point 
vertex function: 

e_ [{(^e}]^Aa{l+l5he2+pe2ei)+Ba{l+l3hei)5e_„e_, 

(D5") 

The same method can be used to constrain the form 
of r2^'^''[{t/)e}]. We now derive Eq. (|6ip with respect to 
and $6(2/, ^2) and evaluate the result in a field 
configuration $e = 0e- The part that comes with no 
kronecker 5ab satisfies 



do 



(i-/?^i0i)r 



(2) 



0. 



(D6) 



By using the symmetry under permutation of the 2- 

— (2) 

points vertex function, we can parametrize ^ in the 
following manner: 

+ Cabdld2 + Cbad2dl + -C'(jf,6'i6'i02^2 

(D7) 

where we have again dropped the explicit dependence 
on the Euclidean coordinates. After inserting this 
parametrization in Eq. (jD6[) . we end up with a set of con- 
straints that leads to the following form for the 2-point 
vertex: 



with Aab = Aba- 



2) (D8) 



2. Contribution of the "ultralocal" parts of the 
effective action to the 2-point vertices 

A generic "ultralocal" term with one copy can be writ- 
ten as 

E / / l[{d''x,<^a{x^,e))R{{x,}) (D9) 

where the function R is symmetric under permutation 
of any pair of arguments. The second derivative of the 
previous term with respect to $afei) and $6(212); evalu- 
ated for $e = 0e, Ve G {1, n}, yields a contribution to 
ril,flj{0e}] of the form 

r P 

p{p - + Phei)5e_,,e_, / \{[d''x,<j)a{x^))R[{x^}) 

1=3 

(DIO) 



which contributes to Ba in Eq. (jPSp . 

Consider now a generic 2-copy "ultralocal" term. 



IE 



Y[{d''x,^a{x,,9^))^, 



ab "Si "^2 i=l 



(Dll) 



Y[id%<t>b{y,,e2))s{{x^},{yA) 



in the effective action, with S a symmetric function un- 
der permutation of the arguments in curly brackets. Af- 
ter deriving with respect to ^a{xi,0_i) and ^bivijia) 
and evaluating this derivative in a field configuration 
'I' P. = (/'ejVe G we obtain a contribution to 



^(2) 
{a,e,).(,bM2) 



of the form: 



qrii + pe,di)ii + pe2e2)x 



l[id''x,Mx^))l[id%MyJ))sax^}AyA) 

i=2 j=2 

(D12) 



which contributes to Aab in Eq. (jDSj). 

It is easy to show that the generic "ultralocal" terms 
with two or more copies lead to the same Grassmann 
structure for F'^^ as that found in Eq. (|D10I) . For in- 
stance, the generic 2-copy "ultralocal" term (jPlip also 
contributes to F^^^ with the a Grassmann structure pro- 
portional to {1 + I39i6i)6e-^,g^. Similarly, generic "ultralo- 
cal" terms with three or more copies lead to the same 
Grassmann structure for F*^^) as that found in Eq. (ID12p . 



3. Contribution of "non-ultralocal" parts of the 
effective action to the 2-point vertices 

We now add derivatives with respect to the Grassmann 
variables in the effective action. No more than two fields 
can come with Grassmannian derivatives, otherwise the 
derivative would vanish upon evaluating in a uniform 
field configuration $e = (j)e,'^e- € {l;---,?^} and the term 
would not contribute to the 2-point vertex function (in 
such configurations). Without loss of generality, we only 
need to consider the three following 1-copy terms in the 
effective action: 

E/ / (J[d''x,Mx^,0)j J d^yAg_^a{y,0)Tiy,{x.}), 
Y.jJ \ J{d''x,'^>a{x^,9) \ j d%d%il+l3e9)x 



^e'^?a{yl,o)^e^aiy2,0)ui{yuy2}Ax^}) 



and 

E 



J J l'Y[d''x,'fa{x,,e)] J dSid%Ae_^a{yi,0)x 



Ae_^a{y2,0)vi{yuy2}Ax^})), 
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where T, U, V are symmetric functions under permuta- 
tion of the arguments in curly brackets. When deriving 
with respect to 4>a(yi,^i) and $b(2/2,£2) and evaluating 
this derivative in a field configuration {$e = (p^}, these 

terms lead to the following contributions to ^ g [{(f)^}]: 



and 



T(yi,{y2,a;J), 
2(1 + ^9102 + 13929^) 

U{{yi,y2},{xt}) 



and 



SPil + I3h92 + PhOi 
^({yi,2/2},{a;J), 



\{d''x,(t)a{x^)\x (D13) 
\i=2 I 



X (D14) 

^4=1 / 



y fj]d'^a;,</.a(a;Oj x (D15) 



all of them coming with the same Grassmann structure 

(l+^^102+/3^20l). 

— (2) 

Finally we study the contributions to F^^ ^ ) (b e ) [{'^e}] 
of the "non-ultralocal" terms in the effective action. Only 
terms with Grassmannian derivatives on at most two 
fields can contribute. Moreover, in the case where two 
fields arc derived, these derivatives must act on fields in 
different copies. Consequently, we only need to consider 
two kinds of terms: 



ah •'^1 •'I2 •' / yj = i 

(D16) 

/ d''zAg^M^,e,)W{{x,},{y,},z) 



ah "'^1 "'^2 / \j=l 

J d'*zid''z2Ae^$,(zi,^i)Ae^$,,(z2,^2)x 

X({x,},{y,},zi,Z2). (D17) 



After deriving with respect to ^a{xi,9i) and $6(2/1,^2) 
and evaluating in a uniform configuration $e = 0e7Ve € 
{1, n}, one can readily check that these terms give no 
contribution to F^^^. 

This concludes our proof that: 



"ultralocal" parts of the effective action contribute 



toF 



(2) 



ae e by terms proportional to {l+P9i9i)Sg^^g^ 
and to f j^Jg^) (^^^ 



by terms proportional to (1 + 

/3^10l)(l+/3^202). 

• "non-ultralocal" parts of the effective action con- 



tribute to F 



(2) 

a,e, e. 



by terms proportional to (1 



(2) 



P9i92+/3929i) and do not contribute to Fr 
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